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Tables of Certain Functions Occurring in 
Dynamics of Structures 


The following is a short description of tables contained in K. HonEn- 
EMSER and W. PRAGER, Dynamik der Stabwerke. Eine Schwingungslehre 
fiir Bauingenieure Berlin, Springer, 1933.1 vi, 368 p. 15.5 X 23.1 cm. 
These tables were computed with the view to facilitating the solution of 
certain standard problems of Dynamics of Structures. 

1. The frequency functions. The natural frequencies of the small flexural 
vibrations of a prismatic elastic beam which is elastically supported and 
elastically clamped at both ends are found from the roots A of the following 
equation: ? 


eofolesfiD(A) — e:A(A) + fiB(A) + E(A)] 
+ — exS(A) — + A(A)] 
— folesf:B(A) + e:C(A) — frS(A) + 
+ efi:E(A) — 4:A(A) + fiBQ) + DA) = 0, (1) 


A(A) = cosh \ sin \ + sinh A cos A, 

B(A) = cosh \ sin A — sinh \ cos A, 

C(A) = 2 cosh X sin X, 

S(A) = 2 sinh sin 

D(a) = cosh A sin A — 1, 

E(A) = cosh \ cos \ + 1. (2) 


The dimensionless parameters éo, ¢; characterize the rigidity of the clamping 
and fo, f: the rigidity of the supports, the subscripts » and ; referring to the 
two ends of the beam. The roots A of the equation (1) are connected with 
the (circular) frequency w of the flexural vibrations by A = lw'/*(u/ED)™*, 
where / is the length, » the mass per unit length and EI the flexural stiffness 
of the beam. Since any member of a plane frame consisting of straight 
beams may be considered as a beam elastically supported and elastically 
clamped at its two ends, equation (1) can be applied to find the natural fre- 
quencies of any such frame. Experience shows that in such frequency 
computations not only the functions (2) but also the ratios S(A)/B(A), 
C(A)/A(A), and B(A)/D(A) occur frequently. The functions (2) as well as 
these ratios are tabulated, along with 1/A and cotA, p. 323-343, for » 
= [0.00(0.02)10.00; 5D]. These tables were computed by F. W. WALTKING 
and Mrs. K. HOHENEMSER. 

2. The amplitudes of displacement and bending moment in the forced 
vibrations of a beam on two supports. Consider a prismatic elastic beam of 
length / extending along the x-axis of a system of rectangular coordinates 
x, y from x = 0 to x = / and resting on pin-supports at these two points. 
If at the end x = 0 this beam is loaded by a periodic couple C = Cy cos wt 
of a given (circular) frequency the bending moment at the cross section x 
can be computed according to 


t) Cox(a, cos wt, 
101 


102 TABLES OF CERTAIN FUNCTIONS IN DYNAMICS OF STRUCTURES 
where £ = x/l, while \ is defined in the same way as above and 


sinh — . sind(1 — 4) 


The deflection at the cross section x is given by 


where 
_ sinhaA(l — £) , sind(i — £) 
20,9 = 2 sinh 2sink (4) 


The functions (3) and (4) are tabulated, to 4D (p. 344-353) for \ = 0.0(0.1)- 
10.0, and ¢ = 0.0(0.1)1.0. These tables were computed by W. PRAGER. 

3. Tables for the coefficients appearing in the four moment equation of 
structural dynamics. An important problem of structural dynamics is to 
determine the displacements, shear forces and bending moments of the 
members of a plane frame which is loaded by periodic forces of the type 
F = Fycoswt. Once the displacements and bending moments at the two 
ends of a member are found the displacement and bending moments at other 
points of this member can be computed by means of the tables discussed 
in section 2. For the determination of the displacements and bending 
moments at the endpoints of all members of the frame equations have been 
developed * which correspond entirely to the four moment equations of 
structural statics. Whereas in the case of statics coefficients appearing in 
these equations depend only on the dimensions of the structure and on the 
loads applied to it, in the case of dynamics these coefficients depend also 
on the frequency of the periodic loads. It can be shown that all these 
coefficients may be expressed in terms of the following functions: 


#0) = (tanh — tan}), #0) = (tanh + tan}), 


= (coth — cot A), = * (coth cot A), 


= — (cach — ese), VOX) = (esch + ese). 


These functions are tabulated, along with tan (A/2), sec (A/2), cot A, and 
esc A, for 4 = [0.00(0.05)10.00; 5D]. The tables were computed by S. 
GRADSTEIN. 
W. PRAGER 
Brown University 


1Of this work there was a Ukrainian translation, Dinamika Sporudzheni2, Kharkov- 
Kiev, Naukovo-Tekhnichne Vidavnitstvo Ukraini, 1935, 381 p. 
W. P: rager, “‘Zur Berechnung der ebener und raumlicher 
Stabwerke,” Z. te techn. Physik, v. 10, 1929, p. 275-280; F bong ‘Zur Ermittlung 
der Eigenschwingu ngszahlen ebener’ Stabwerke, - Ingenieur-Archiv, v v. 2, 1931, p. 247-274. 
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F. Reinitzhuber, “Beitrag zur Berechnung der Eigenschwingungszahlen raumlicher Stab- 
werke,” idem, v. 8, 1937, p. 349-363. 

3 rager, “Die Beans pruchung von Tragwerken durch schwingenden Lasten,” 
Ingenieur-Archiv, v. 1, 1930, p. 527-532; S. Gradstein and W. Prager, * ‘Beanspruchung und 
Formanderung von Stabwerken bei erzwungenen Schwingungen,” Ingenieur-Archiv, v. 2, 
1931, p. 622-650; “‘Berichtigungen,” v. 3, 1932, p. 434. Tables of tan 4/2, cot A, csc i, 
sec \/2, are given on p. 643-646; and those for the last group of functions listed above, 
p. 64 647-650. See also E. Fliegel, ‘‘Die Elasticitatsgleichungen zweiter Art der Stabwerks- 
dynamik,” idem, v. 9, 1938, p. 20-28, where ¢/(¢? + ¥), we x ¥), 1/(@ — ¥), and 1/¢ 
are tabulated for \ = [0.00(0.50) 5. 00; 2D]; and R. E. Gaskell, ‘On moment balancing in 
mg eel Quarterly of Applied Mathematics, =. University, v. 1, 1943, 
p. 237-249. 
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108[K].—W. G. Cocuran, “The distribution of the largest of a set of esti- 
mated variances as a fraction of their total,” Annals of Eugenics, v. 11, 
1941, p. 47-52. 20.2 X 27.5 cm. 


In this paper the distribution of 
g = (largest 


is studied for the case where 1, “2, +++, 4, ***, %, are mutually independent random vari- 
ables each distributed as x?0* for r degrees of freedom, i.e., 
b(us)dus (204 du, (¢=1,2,---,m). 
The values of g which have a probability 0.05 of being exceeded in random sampling are 
tabulated for » = [3(1)10; 4D], r = 1(1)6(2)10. 

The utility of these “5% points of g” derives from the fact that u may be interpreted 
as the “sum of squares with r degrees of freedom” corresponding to an unbiased estimate, 
s? = u/r, of the variance, o*, of a Gaussian population. In consequence, the distribution of 
g provides the distribution stated in the title, and the “5% points of g” provide a test of 
whether a variance estimate which appears to be anomalously larger than the others of a 
set is an estimate of the same o*. 

The distribution of g for the case of r = 2 was obtained by R. A. Fisher, “‘Tests of 
significance in harmonic analysis,” R. So. London, Proc. v. 125A, 1929, p. 180-190, who 
proposed its use in the harmonic analysis of a series to test the statistical significance of a 
particular term which is picked out on inspection because of its exceptional magnitude. In 
this paper Fisher gives a table of the 5% points of g to 5D for r = 2 and m = 5(5)50. In 
R. A. Fisher, ‘The sampling distribution of some statistics obtained from non-linear equa- 
tions,” Annals of Eugenics, v. 9, 1939, p. 238-249, the distribution of g is employed to illus- 
trate the disturbance introduced into the analysis of variance test of a regression equation 
when the equation is ‘non-linear’ in the sense that the eliminant of the regression coefficients 
is a non-linear function of the observations which satisfy the regression equation exactly. 
W. L. Stevens, “Solution to a geometrical problem in probability,” Annals of Eugenics, v. 9, 
1939, p. 315-320, found that the distribution of g for r = 2 provides the answer to the 
following question: ‘‘On the circumference of a circle of unit length, ” arcs, each of length x, 
are marked off at random. What is the probability that every point of the circle is included 
in at least one of the arcs?” This curious coincidence is explained by R. A. Fisher, “On the 
similarity of the distributions found for the test of significance in harmonic analysis, and 
in Stevens’ problem in geometrical probability,” Annals of Eugenics, v. 10, 1940, p. 14-17, 
and the distribution of the second largest ratio, gz, when r = 2 is also discussed. This paper 
is followed, p. 17, by W. L. Stevens, “‘A table of the test of significance in harmonic analysis,” 
in which the 5% points of g: (our g) and gz are tabulated forr = 2and m = [3(1)10(5)50; SD]. 
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The only full table of the distribution of g which I have come across is Table 2 in H. T. 
Davis, The Analysis of Economic Time Series, Bloomington, Indiana, Principia Press, 1941, 
where for r = 2 the probability of g exceeding k/n in random sampling is tabulated to 5D 
for k = 0.1(0.1)10.0; m =°10(10)70 and for k = 5.1(0.1)10.0; m = 80(10)160(20)300. 

C. EIsENHART 


109[K].—F riepa S. Swep and C. EIsEnnart, “Tables for testing random- 
ness of grouping in a sequence of alternations,” Annals Math. Stat., 
v. 14, 1943, p. 66-87. 16.3 X 23.4 cm. 


If two different kinds of objects are arranged along a line, they will form two or more 
distinct sets of like objects. For example, if 4 a’s and 3 b’s are arranged as follows: aaabbab 
there are 4 such groups. In general, if there are m objects of one kind and n of another, there 
are C.*" different possible arrangements. The proportion of these arrangements having u 
distinct sets of like objects is 


(1) fu = , if w=2k, even) 
cat 
Ga + if u=2k—-1, odd). 


, 


We can assume m < n without loss of generality. Then (1) holds for k = 1, 2, ---,m+ 1. 
The probability law (1) was derived by W. L. Stevens (Annals of Eugenics, v. 9, 1939, 
p. 10-17). Wald and Wolfowitz (Annals Math. Stat., v. 11, 1940, p. 147-162) have used 
the function x as a criterion for testing the statistical hypothesis that a sample of m elements 
and a second independently drawn sample of elements (both known to be from continuous 
probability distributions) could have come from identical probability distributions. Wald 
and Wolfowitz have shown that if m‘and n are large, then the quantity 


2 )/2(<= 
where a = m/n, is approximately distributed according to the normal probability law 
1 
— 


Now let us write 
Plu = 


The authors have presented in Table I the values of P{u = u’'}, to 7D, for m=n = 20 
and 2=m= 20. 

Now let us denote by 1, the largest integer for which P{u = u,.} =e when e < 0.50, 
and the smallest integer for which P{u < u.} => when e > 0.50. Table II gives the values 
of u, for « = 0.005, 0.01, 0.025, 0.05, 0.95, 0.975, 0.99, and 0.995, and for m = n = 20 and 
2=m= 20. 

By making use of the approximately normally distributed variable x in (2) the authors 
have given in Table III values of u, for m = n from 10 to 100. In making this tabulation, 
a continuity correction was used in order to improve the approximation to u, for small 
values of m and n; in particular the values 10 <= m = n < 20. The authors give several 
examples to illustrate the use of their tables. 

These tables should prove extremely useful as a basis for making simple significance 
tests in a wide variety of statistical problems which can be reduced to a treatment of 
sequences of two kinds of elements. 

S. S. W. 
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110[K].—E. S. Pearson and H. O. Hart ey, “The probability integral of 
the range in samples of 1 observations from a normal population,” 
Biometrika, v. 32, 1942, p. 301-310. 19.3 K 27.3 cm. 


If x1, x2, ---, X, denote a random sample of m observations arranged in ascending order 

of magnitude, the sample range is defined as x, — x1. The authors consider the quantity 
Xn — 
y= 
assuming that the sample has been drawn from a normal or Gaussian population with 
standard deviation c. If f,(w)dw is the probability element of w, the cumulative distribution 
law of w is defined as 
PAW) = J.” 

It is shown that 


PAW) = + am f™ ( a, 


a(x) = 


In Table 1 of the present paper, values of P,(W) are calculated, to 4D, for W = 0.00 
(0.05)7.25 and for m = 2(1)20. 


L. H. C. Tippett (Biometrika, v. 17, 1925, p. 386-387) has tabulated the mean value of 
w (ic. the value of f° whalw)dw ) for m = 2(1)1000. 


If a constant a, is chosen so that the mean value of an(x, — x1) is equal to o, then 
@,(X%n — 1) may be considered as an estimate of ¢, although not a highly efficient one. The 
efficiency of this estimate of ¢ has been discussed by O. L. Davies and E. S. Pearson (R. 
Stat. So., Jn., supp. 1, 1934, p. 76). In Table 2 of the paper under review, values of a, are 
tabulated for = 2(1)12. In Table 2 are also tabulated values of W such that P,(W) = 0.001, 
0.005, 0.01, 0.025, 0.05, 0.10, 0.90, 0.95, 0.975, 0.99, 0.995, 0.999 for m = 2(1)12. 

S. S. W. 


where 


111[K].—MaxinE MErRINGTON, “Table of percentage points of the 
t-distribution,” Biometrika, v. 32, 1942, p. 311. 19.3 & 27.3 cm. 


By setting m; = 1 and m2 = n in the incomplete beta-function 
+ 
t= — x)/x 
one obtains the distribution function 
2r(3n + 3) dt o< 
Vn Vx (1 + 


of the absolute value of the ‘‘Student” ratio for m degrees of freedom. Mrs. Merrington has 
tabulated values of tp for which 
Si = P, 


to 5D, for P = 0.50 ,0.25, 0.10, 0.05, 0.025, 0.01, 0.005, and for m = 1(1)30, 40, 60, 120, ~. 
The table was derived from Miss Thompson's table RMT 99. 


— x) 


and setting 


gn(t)dt t< 


S. S. W. 
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112[K].—E. S. Pearson and H. O. Hart ey, “Tables of the probability 
integral of the studentized range,” Biometrika, v. 33, 1943, p. 89-99. 
19.3 X 27.3 cm. 


Denote by x1, x2, --*, X, a random sample of m observations (in ascending order of 
magnitude) from a normal population with standard deviation c. The authors of the present 
paper have tabulated (RMT 110) the probability integral of the quantity w = (x, — x1)/c. 

In most practical statistical problems ¢ is unknown and hence the ratio w is unusable. 
Hence, it is desirable in problems of this kind to replace « by an estimate of ¢ calculated 
from a sample. In the present paper, the authors consider two independent samples, the 
first consisting of m observations from which x, — x; is calculated, and the second consisting 
of » + 1 observations, X1, X2, ---, X»4: from which an estimate s of o is calculated as 
foilows: 

1 
(1) 2(X;— X), 
where X is the arithmetic mean of the » + 1 observations of the second sample. Instead of 
w the quantity 
Xn — X11 
s 


is considered. Denoting the probability law of g by gn(q)dg, and the probability integral by 


(2) »P 'n(Q) v8n(q)dq, 
an approximation of the following form for »P,(Q) is found 
Pa(Q) = Pa(Q) + ~a0(Q) + = ba(Q) 
where P,(Q) is the probability integral of w and 
1 
1 (Q'd'P, QP, QP, 
= 2 aw saw 2am Paw 


P,(Q), an(Q), and b,(Q) are tabulated to four, two and one places of decimals respectively 
for Q = .00(.25)6.00, » = 3(1)20. The authors do not guarantee the last figure in 6,(Q) for 
values of m in the upper end of the range. Values of Q are tabulated to two places of decimals 
(in some cases only one) for »P»(Q) = 0.01, 0.05, 0.95, 0.99 using approximation (3), for 


n = 2(1)20, vy = 10(1)20, 24, 30, 40, 60, 120, ~. 
Harmonic interpolation with respect to » is made possible by the tabulations for this 


particular set of values of ». 
S. S. W. 


113[B, D, E, G, H, L]|—Paut RupoL_F EuGEN JAHNKE (1863-1921) and 
Fritz Empe (1873- ), Tables of Functions. (|Funktionentafeln mit 
Formeln und Kurven| . . . Published and distributed in the public interest 
by authority of the Alien Property Custodian under license No. A.57. 
New York, Dover Publications, 31 East 27th Street, 1943. [xv], 303, 
76 p. 15.4 XK 29.8cm. $3.50. 

There is hardly any single volume which could be more useful to the general applied 


mathematician and the mathematical physicist, than this remarkable work of Jahnke and 
Emde, which has been developed during the past 34 years. The Alien Property Custodian 
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is to be congratulated not only upon his choice of such a work for reproduction at the present 
time, but also upon the publication price which places the volume within the reach of ordi- 
nary individuals and small libraries, as well as of those more opulent. That the opportunity 
was appreciated is illustrated by the fact that a thousand copies were ordered before 
publication on June 1. 

The work started out in modest fashion as one of the series Mathematisch-Physikalische 
Schriften fiir Ingenieure und Studierende, no. 5, ed. by E. JAHNKE, and with the title 
1. Funktionentafeln mit Formeln und Kurven, Leipzig, 1909, xii, 176 p. 15.3 X 23.7 cm. With 
53 text figures. The text contained the following XIII main headings: I, The functions 
x tanx and tanx/x; II, Roots of transcendental equations; III, Transformation of a + bi 
into re**, and conversely; IV, The exponential functions e* and e~*; V, Hyperbolic functions; 
VI, Sine, cosine, and logarithmic integral; VII, Fresnel integrals; VIII, Gamma function; 
IX, The Gauss error integral ¢x; X, The Pearson function F(r, v); XI, Elliptic integrals _ 
and functions (p. 46-79); XII, Spherical harmonics (p. 79-89); XIII, Bessel functions (p. 
90-174, that is, about half of the volume). The work is filled with formulae, and numerical 
tables and a notable collection of more than fifty interesting text-figures, from which not 
only the general variation of the functions may be observed, but also their approximate 
numerical values can often be readily estimated for given arguments. There are many 
references to the literature, where more extensive material is to be found. The tables given 
are usually short and to a small number of decimal places. When this volume was published 
no tables existed for many of the higher functions, but after the first world war great ad- 
vances had been made in this regard, so that there was an embarrassment of riches to deal 
with when EmpE (his colleague having died) came to prepare the second edition nearly 25 
years after the first one had been published. This volume, nearly double that of the first 
edition in size, was in two languages throughout—German and English, with two title pages, 


2. Tables of Functions with Formulae and Curves, second (revised) edition, with 171 figures. 
Funktionentafeln mit Formeln und Kurven, Leipzig, 1933. xviii, 330 p. 16 X 24.2 cm. Compare 
RMT 32. 

The number of sections is here increased from XIII to XIX, ten of which were devoted 
to the elementary functions, the number of pages devoted to them being increased from 18 
to 78. These sections bore the following titles: I, Table of powers; II, Auxiliary tables for 
computation with complex numbers; III, Cubic equations (graphical and tabular aids to 
their solution); IV, Elementary transcendental equations; V, x tan x, tan x/x, and sin x/x; 
VI, Exponential function; VII, Planck’s radiation function; VIII, Source functions of heat 
conduction; IX, Hyperbolic functions; X, Circular and hyperbolic functions of a com- 
plex variable. These sections were concluded with an Index of tables of the elementary 
transcendentals. 

Then comes XI, Sine, cosine, and logarithmic integral tables. The Pearson function 
table is eliminated (see MTAC, p. 70), and the Gauss error integral and Fresnel integrals 
are brought together under XIII, Error integral and related functions. The Gamma func- 
tion is here treated under XII, Factorial function, much enlarged and greatly enriched with 
charts and relief drawings. Section XIV, Theta functions, is now separated from the single 
section on elliptic integrals and functions which has been divided into XV, Elliptic integrals, 
and XVI, Elliptic functions. The material here has been much expanded and the drawings 
and reliefs are indeed notable new features. So also with respect to XVII, Legendre func- 
tions, and XVIII, Bessel functions (p. 192-318). The latter section has been greatly ex- 
tended. The Nielsen definition of the functions of the second and third kind (N and H— 
although Y and H are more commonly used in Great Britain, but this Y is not used by 
Watson) is now uniformly introduced in all of the tables. The final section, XIX, The Rie- 
mann zeta-function, is new. In every section after X the lists of “text-books and other 
tables,” adds largely to the usefulness and value of the volume. The representation of com- 
plex functions by surfaces whose ordinates are equal to the modulus of the function, have 
been referred to above as “reliefs” of the function. This graphical method, adopted here 
for the first time, is an innovation of outstanding importance. 
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Even without the introductory article of our last issue special reference would naturally 
be made to the following quotation from the preface: “As in the first edition great use has 
been made of the work of the Mathematical Tables Committee of the British Association. 
Fortunately this committee has decided to publish collections of the very accurate tables 
which they have calculated ia past years. Two volumes have already been published.! The 
mathematicians, physicists, and engineers of the whole world regard with the greatest 
wonderment and gratitude this colossal undertaking of their English colleagues, who have 
taken upon themselves almost entirely the heavy load of new computation. It is hardly to 
be conceived that other countries can continue much longer to look idly on without helping 
in this work.” 

Empe published 3. third (revised) edition of his work, “with 181 figures,” in Leipzig, 
Teubner, 1938, xii, 305 p. Here again important new revisions and additions were made, 
especially with a view to providing advanced research workers with an up-to-date desirable 
tool. Emde summarizes as follows the chief points in which the third edition differs from the 
second: “In the complete elliptic integrals of the first and second kind, formulae and numeri- 
cal tables have been added for other than the Legendre standard forms. The numerical calcu- 
lation is thereby improved in many cases. In the cylinder functions the Debeye series have 
been brought into a more convenient form by limiting them to real values of the index and 
by limiting the angle of the argument to values which occur in practice; they have thus been 
freed from the additional determinations which necessarily added to the difficulties of the 
general case. Detailed numerical tables are given for the Lommel-Weber and Struve func- 
tions of orders zero and one. Formulae and graphical representations are given for the con- 
fluent hypergeometric functions and for the Mathieu functions associated with the elliptic 
cylinder, which will enable these two classes of functions to be employed more widely in 
scientific calculations.” 

This new edition had also 40 new figures and reliefs. The elementary functions occupying 
76 pages of the second edition are not included in the third.2? These pages have been 
elaborated into an independent work, 


3A. oF. Empe, Tafeln Elementarer Funktionen, Leipzig, Teubner, 1941, xii, 181 p., 83 figures, 
text in German and English. I know of no copy of this work which has reached this country. 
The eleven main topics of no. 3 are now as follows: I, Sine, cosine, and logarithmic integral; 
II, Factorial function; III, Error integral and related functions; IV, Theta-functions; V, 
Elliptic integrals; VI, Elliptic functions; VII, Legendre functions; VIII, Bessel functions 
(p. 126-168); IX, The Riemann zeta function; X, Confluent hypergeometric functions; 
XI, Mathieu functions. 

4. Of this third edition in 1941 G. E. Stechert and Co. of New York, purchased from 
Teubner the right to print a facsimile edition. Since this reprint (now sold at $3.00) has been 
constantly available for purchase by scholars since 1941 it caused considerable surprise to 
learn that in 1943 the Alien Property Custodian had authorized the Dover Publications 
to produce 
5. the fifth edition of JAHNKE and EmpeE’s work, the volume under review. This volume 
includes not only no. 3, but also the first 76 pages (about elementary functions) of no. 2. 
In most respects the volume makes an excellent impression. 

We have good reason for believing that while in ignorance of the Stechert edition the 
Alien Property Custodian was fully justified in authorizing the production of no. 5. As we 
have indicated above, no. 4 is also for sale, until the edition is exhausted. 

References to errors may now be mentioned. Errors in F(¢, a) include, (a) F(60°, 60°), 
for 1.2125 read 1.2126; (b) F(37°, 75°) for 0.6913 read 0.6919. These are in no. 1, p. 57, 58; 
in no. 2, p. 137, 138; in nos. 3-5, p. 65, 66. Compare RMT 50 (Scripta Mathematica, v. 3, 
p. 365). Mr. W. D. LamBeErt, of the Coast and Geodetic Survey, discovered an error 


in no. 2, p. 131, and in nos. 3-5, p. 59, in one of three forms of elliptic integrals equal 


dt 
to 4F(k, ¢). “In the first of these fi = 
4F(k, ¢). “In the ese forms, cos 1 + x*. The corre- 
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sponding modulus should then be 1/72 or sin 45°, and not 2(¥2 — 1)¥2 = 0.985 1714 
dt 
= sin 80°.1207. This latter value of the modulus, ;, is obtained f {az- 
us, Ri, 1S al rom 
(2 — V2) F(w, k:), where w is determined from 
tanw = + 1)(1 — x)/(1 +2). 


These are the familiar Landen relations k; = 2Vk/(1 + &) and sin (20 — ¢) = k sing.” 

Mr. LAMBERT pointed out to me also that O. S. Apams in his Elliptic Functions Applied 
to Conformal World Maps, Washington, 1925, p. 58, recorded two serious errors in the table 
(to 4D) of Weierstrass functions, no.1, p.75, namely: the value of g(u) for r = 3 should be 
1537.9625 instead of 1468.820; and the value of g’(u) for r = 35 should be — 75.9603 in- 
stead of — 73.4302. These errors occur also in all later editions of J. and E. as weil as in the 
table computed by A. G. Hapcock (to 5D) and published by A. G. GREENHILL in R. Artil- 
lery Institution, Minutes and Proc., v. 17, 1889, p. 212-213. Mr. Adams verified the accuracy 
of other values of g(u) and g’(u) in this table for r = 5 to r = 59 inclusive; it continues 
tor = 240. 

In Amer. Math. So., Bull., v. 3, 1897, p. 153-155, R. W. WILLson and B. O. Perrce 
give a “Table of the first forty roots of the Bessel equation Jo(x) = 0 [to 10D] with corre- 
sponding values of J:(x) [to 8D].”” J:(17) is given as +0.10999114;. In all editions of J. and 
E. this is abridged to 4D. Harry BATEMAN pointed out to me that for J:(17), +0.1010 is 
always given; instead of +0.1100. He drew my attention also to the fact that in the table of 
the roots of Jo(x) No(kx) — Jo(kx)No(x) = 0, on p. 204 of nos. 3-5, p. 274 of no. 2, and p. 162 
of no. 1, all 6 roots for each of three values of k (1.2, 1.5, 2.0) are wrong. These 18 values 
are given correctly in M. Musxat, F. Morcan, M. W. Meress, “The lubrication of plane 
sliders of finite width,” Jn. Appl. Physics, v. 11, 1940, p. 212. The table printed in J. and E. 
first appeared in A. KALAHNE, Zeits. f. Math. u. Phys., v. 54, p. 68, 1906. 

Mr. R. D. Brown Jr., of Blue Bell, Penna., reported that the value of J,(21) = — 0.3175, 
given on p. 248 of no. 2, and p. 177 of nos. 3-5, does not agree with what appears to be the 
correct value, —.03175 (verified by H. Bateman), given on p. 296 of A. Gray and G. B. 
Matuews, A Treatise on Bessel Functions, second ed. by A. Gray and T. M. Macrosert, 
London, 1922. 

We shall be happy to receive for publication other reports of errors in this wih. 


1 Emde’s reference is doubtless to B.A.A.S., Mathematical Tables, v. 1, Circular and 
Hyperbolic Functions, Exponential Sine and Cosine Integrals, Factorial (Gamma) and Derived 
Functions, Integrals of Probability Integral, London, 1931; and v. 2, Emden Functions being 
Solutions of Emden’s Equation together with Certain Associated Functions, London, 1932. 
The first volume published by the Committee was, however, A. J. C. Cunningham’ sA 
Binary Canon, London, 1900, viii, 172 p. In 1901 the Committee recommended the publi- 


cation of the same author's Quadratic Partitions, but this work was published otherwise in 
1904.—EpirTor. 


2 The “Index of tables of the elementary transcendentals,” occupying P. 76-78 of = 
second edition, has been revised and improved on p. 299-301 of the third edition. 

date of publication (1921) of U. Meyer, Fluchtlinientafeln des Hyperbeltangens einer eee 
Verdnderlichen is, however, still lacking. 


114[Q, U].—A. A. Aceton, Manual of Celestial Navigation, New York, Van 
Nostrand, 1943. 103 p. 17.4 X 244 cm. $1.85. The first printing 
of this book was in April 1942; the fifth in February 1943. 


This volume contains two basic tables which allow one to use an Assumed Position 
method, similar to but an improvement over that of J. Y. Dreisonstok, Navigation Tables 
for Mariners and Aviators, (H.O. 208, RMT 103), or a Dead Reckoning Position method 
identical with that presented by the author in his Dead Reckoning Altitude and Azimuth 
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Table (H. O. 211, RMT 104). Both methods give an accuracy of solution to within five- 
tenths of a minute of arc in altitude without interpolation; this is adequate for all ordinary 
navigation on land, at sea and in the air. The author is a Commander in the U. S. Navy 
and Head of the Navigation Division of the Department of Seamanship and Navigation at 
the United States Naval Academy at Annapolis. 

In the notation of RMT 103, in his AP method, Ageton drops a perpendicular from Z 
on the great circle through P and S. The foot of the perpendicular is X and the sides PX 
and ZX are called 90° — K and R respectively. Table I is a double-entry table in which are 
tabulated three quantities associated with the polar triangle, B, K, and Z’, for integral 
values of the arguments L(0° to 90°) and #(1° to 179°). K has already been defined as the 
declination of X; it is given to the nearest tenth of a minute of arc. B is 10° log csc R, and 
is given to the nearest integer except that one decimal is given for L(0° to 63°) and #(1° to 
20°), and for L(64° to 90°) and #(1° to 16°). Z’ is the angle PZX and is given to the nearest 
tenth of a degree. All of the values corresponding to a single latitude are presented in two 
columns on two pages facing each other, and may be reached with a single opening of the 
book! This marks a great improvement over H. O. 208. 

The formulas used to complete the solution are 


A(t) = + B(K~d), Z2Z=2'+2". 


As in H. O. 211, (K ~ d) is used to indicate] K — d|, and A(x) and B(x) are used to denote 
10° log csc x and 105 log sec x respectively. 

Table II is similar to the one contained in H. O. 211 except that only the values of A 
and B corresponding to angles with an integral number of minutes of arc are given in their 
entirety; the end figures which differ from the complete values above them are given opposite 
the half-minutes. This makes the table seem somewhat more open, and perhaps easier to 
read. It should be mentioned that the printing of the second table is poor; this is especially 
noticeable if one has been using the Hydrographic Office Tables which are printed at the 
Government Printing Office. 

Both Tables I and II have good thumb indices along the right hand edge of the page. 
They are easy to use but would be even easier to use if the back cover of the book were 
more flexible. In the index to Table II, the value of A given beneath 0° is 38400 and should 
be 384000. This will doubtless be corrected in another printing. 

Inasmuch as the DRP method given in this volume is identical with that presented in 
H. O. 211 (RMT 104), there seems to be no point in repeating the discussion of its ad- 
vantages and disadvantages. 

A comparison of Table I with a corrected copy of Dreisonstok, Navigation Tables for 
Mariners and Aviators (H. O. 208) indicates that the Manual has a much smaller number of 
large errors. Some of its errors are obviously inherited. Ageton mentions in the Preface that 
he has studied the works of Souillagouet, Bertin, Ogura, Weems, Dreisonstok, and Gingrich 
before undertaking the construction and arrangement of Table I. No mention is made of 
Hughes’ Tables for Sea and Air Navigation by L. J. Comrie which first appeared in 1938 
and which contains material similar in content and arrangement to that in the Manual; 
compare RMT 115. Ageton’s method of determining azimuth is different from that of Comrie. 

It is interesting to note that the beautifully simple and rapid short method of computing 
azimuth included in the first printing is not contained in later printings. In a letter, the 
author states that it ‘‘proved to have unpredictable errors in it and so had to be eliminated 
from the later printings.” 

This is an excellent book to use in the teaching of celestial navigation. The elegance 
and simplicity of its two methods reduce the tedious details of a computation. Since it 
contains the elements for both an AP and a DRP method, it should appeal to many users. 
One suspects that a student well trained in the use of this volume will prefer it to any other. 

C. H. 
Brown University 
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115[Q, U].—L. J. Comrie, Hughes’ Tables for Sea and Air Navigation, 
London, Henry Hughes, 1938. lIvi,182p. 15 X 24.6cm. 20 shillings. 


Three separate tables and two methods, one AP and the other DRP, are offered for the 
solution of the fundamental problem of celestial navigation. An accuracy of approximately 
0’.4 in computed altitude can be attained without interpolation by the use of the AP method; 
by using the DRP method with interpolation, an accuracy of 0’.1 can be had. The former is 
adequate for most navigation, more than adequate for air navigation; the latter will meet 
the most rigid demands of navigators. Comrie was formerly the Superintendent of the 
British Nautical Almanac Office and is one of the outstanding authorities of the world on 
the design, construction and accuracy of tables. 

— Using again the notation of RMT 103, Comrie chooses to divide the astronomical tri- 
angle into two right triangles by dropping a perpendicular from Z on the great circle through 
P and S. Call the foot of this perpendicular, X, (Comrie uses Q), and let PX = 90° — K 
and ZX = R. Also let 

B(x) = 10° log sec x 
C(x) = 10° log csc x 


In Table I are given four quantities associated with the polar right triangles, K, A, D, 
and Z,, with integral degree arguments, L(0° to 89°) and #(0° to 180°). As indicated above, 
K is the declination of the point, X; it is given to the nearest tenth of a minute of arc. A 
is B(R); it is given to the nearest integer except that for values of A less than 665, and L 
between 0° and 34° inclusive, one decimal is given. D is 10-*C(R); it is given to the nearest 
integer. Z, is the angle PZX;; it is given to the nearest tenth of a degree. 

Table I is undoubtedly the most accurate table of this type yet available to the general 
public, although there are approximately 500 insignificant errors of one unit in the last 
place. Comrie recomputed the values for all latitudes 62° to 89°; this section is remarkably 
free from error. Since all of the values corresponding to a single latitude are assembled on a 
single page, the table is extraordinarily convenient to use. No edge-thumb-index of the 
table is offered, although notched pages allow one to turn at once to each of the various 
tables. 

The formulas used in the AP method are: 


Chk) =A +B(K~d), 10*log tan Z; = D + 10*log tan(K ~d); Z= Zi: 


Table II isa table of B(x) and C(x) with argument x = [0°(0’.5)180°]. Band C are given 
to the nearest integer except that one decimal is given for all values less than 665. Comrie 
has carefully eliminated rounding-off errors and the values given are correct to the last 
place. 

Table III is a table of values of 10* log tan x, each given to the nearest integer, with 
argument x = [0°.1(0°.1)89°.9]. The number — 10000 is omitted for values of x less than 
45°. This entire table is printed on two pages facing each other. 

The formulas used in Comrie’s DRP method are: 


C(R) = B(L) + Cd, C(K) = C(L) — B(R), 
C(h) = B(R) + B(K~d), 108 log tan Z; = 10-*C(L) — 10* log tanh, 
10? log tan Zz = 10-*C(R) + 10* log tan (K — d). 


Some persons will object to the elaborate and detailed explanations which occupy 56 
pages at the beginning of the book. Others will appreciate the security afforded by having 
a complete explanation permanently attached to the tables. 

In the Preface, Comrie modestly claims originality only in arrangement, accuracy, 
typography and explanation; it is precisely the high quality of these features that makes 
this book so outstanding in its field. In fact, if one looks for something to criticize, there are 
only these minor defects: three separate tables are required for each of Comrie’s methods; 
no edge-thumb-index is offered for Tables I and II; some shifting of the decimal point is 
required at two places in his DRP method. 


j 
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In addition to the three principal tables, a table of meridional parts to one decimal for 
latitudes [0°(1’)80°] and a traverse table giving Difference in Latitude, and Departure, 
each to one decimal, with arguments Bearing [1°(1°)89°] and Distance [0(1)100(100) 1000) 
are given. Some twelve other auxiliary tables are presented which will meet most of the 
requirements of ordinary navigation. 

C. H. SmiLey 


116[C].—THEopoR WITTSTEIN, Addition and Subtraction Logarithms to Seven 
Decimal Places, Chicago, Ill., Charles T. Powner Co., P.O. Box 796, 
1943. viii, 126 p. 17.2 X 22.9 cm. $2.00. 


The name of the real author of this work is nowhere mentioned in the volume, but 
instead, after the above title, is merely given “Edited by William W. Johnson, Instructor 
in Applied Mathematics John Huntington Polytechnic Institute.” This so-called editor 
contributed nothing but three and one half brief pages of illustrative material, partly 
copied, partly highly misleading, and quite inadequate for instructing the novice in the 
use of the volume. Pages 2—126 are complete facsimile reproductions of exactly these pages 
in the following well-known work of WITTSTEIN, with title pages and introductory matter 
in French and German: Logarithmes de Gauss @ sept Décimales pour servir a trouver le Lo- 
garithme de la Somme ou de la Différence de deux Nombres, leurs Logarithmes étant donnés, 
arrangés d’aprés une nouvelle méthode . . . Ouvrage servant de supplément @ toute table ordi- 
naire de Logarithmes a Sept Décimales. Siebenstellige Gaussische Logarithmen zur Auffindung 
des Logarithmus der Summe oder Differenz zweier Zahlen, deren Logarithmen gegeben sind .. . 
Ein Supplement zu jeder gewéhnlichen Tafel siebenstelliger Logarithmen, Hannover, 1866, 
XVI, 127 p. 17.9 X 25.2 cm. 

On the back of the title page of the work under review is “Copyright, 1943 P.R.C. 
Publications, Inc.”’ For reasons suggested above, it is hardly conceivable that this notice 
implies any real equity on the part of the copyrighter, any more than it did in the case of 
William Charles Mueller who put his own name and the following original title on the cover: 
A New Manual of Natural Trigonometrical Functions to Seven Places of Decimals of Sines and 
Cosines of the Angles from 0 to 10 000 seconds [sic] of a work wholly prepared by W. JoRDAN 
(see MTAC, p. 8), before having it copyrighted in 1907. The title proves that “editor” 
Mueller also did not really understand what was in his volume. America’s record of pirating 
mathematical tables is getting to be unenviable. We learn that the fourth edition of the 
Barlow-Comrie Tables (RMT 82), not copyrighted in this country, has been pirated by the 
Chemical Publishing Co. of Brooklyn, N. Y., some time American agent for the sale of this 
English work. These are but three illustrations. 

Addition and Subtraction Logarithms are used to obtain the logarithm of the sum or 
difference of two numbers when only their logarithms are known. Thus the result is found 
in one table entry, instead of three entries necessary when the ordinary table of logarithms 
is employed. If a and b are two positive numbers, a > 5, then in the notation of the book 
under review 

log a — log b = log (a/b) = B = log x; 
log (a@ — b) = logb +A or A = log (x — 1). 


The table gives values for B or A from which the corresponding A or B are found. Since 
x >1,B>0,and A =0 whenx =2. For this range B lies between 0 and log 2 = 0.301 0300; 
and A lies between — © and 0. But the mantissas of A and B would be the same to 
about 7D at —7; hence the range for A can be taken from —7 to 0, or 3 to 10 if 10 is 
added to these values of A. It is these latter values for A which are given, p. 1-45, where 
A = [3.0(0.1)4.00(0.01)6.000(0.001)8.0000(0.0001) 10.0000; 7D]. 

If log (a/b) = A = log x; and log (a + 5) = log) + B, then B = log (x + 1). Here 
A >0, B > 0.3010300, and when A reaches 7.0, B has become 7.000 0000. Here A = 
[0.0000(0.0001)4.0000(0.01)6.00(0.1)7; 7D], p. 46-126. 
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Wittstein’s original page presented an excellent appearance which is naturally somewhat 
marred in reproduction, slightly reduced in size, with curtailed margins. 

Another seven-place table of Addition and Subtraction Logarithms, by J uLtus Zecu, 
appeared as Table XII (p. 636-836, introduction p. XXI-XXIII) in the 1849 edition of 
Vega’s Sammlung mathematischer Tafeln, edited by J. A. Hiilsse. This was reprinted with 
an independent title page, but its original page numbers, at Leipzig and Berlin in 1849, 1892 
(third) and 1910 (fourth ed.). The Zech form of table is now generally preferred to that of 
Wittstein so that a reproduction of this table would have been decidedly preferable to that 
of Wittstein. A modern 6-place form of the Zech table, is B. Coun, Tables of Addition and 
Subtraction Logarithms with Six Decimals, second ed., London, Scientific Computing Service, 
1939; there is an excellent preface by L. J. Comrie. 

After the above-mentioned 7-place tables of this type, there is a single table with a larger 
number of places, namely: M. H. ANDoyer, “Tables fondamentales pour les logarithmes 
d’addition et de soustraction,”” Bull. Astronomique, s. 2, v.2, 1922, p. 5-32; the final Table III 
(p. 32) is of 10* for x = [0.00(0.01)1.00; 15D]. As in the Zech table, if log a — log b = d, 
log (a + b) = loga + A, log (a@ — b) = loga — S, A = log (1 + 10-4), S = — log (1 — 107%), 
and the table is for the Addition and Subtraction Logarithms A and S in terms of the 
parameter d = [0.00(0.01)9.00; 16D]. 

A truly American table, original in form, is 
L. M. BERKELEy’s Addition-Subtraction Logarithms to Five Decimal Places, New York, White 
Book and Supply Co., 1930. 

The idea of a table of Addition and Subtraction Logarithms was first suggested by 
ZECCHINI LEONELLI (1776-1847) who for a time taught mathematics and architecture at 
Bordeaux. Here in 1802 or 1803 [an onze] he published his 
o Supplément Logarithmique, contenant La décomposition des grandeurs numériques quelconques 
en facteurs finis; . . . et La théorie des logarithmes additionels et déductifs, ou de certains 
logarithmes qui donnent directement les logarithmes des sommes et des différences des valeurs 
naturelles, dont on ne connoit que les logarithmes.' This edition is not listed in the published 
catalogues of the Bibliothéque Nationale, or British Museum, or Library of Congress. 
Leonelli considered that a table of this kind to 14D was desirable. 

Apparently it was this particular edition of Leonelli which first drew Gauss’s attention 
to the subject. In a paper * of 1812 Gauss wrote about the importance of such logarithms 
in practical computing and gave a five-place table where log (1 + 1/x) and log (1 + x) are 
given with their differences for argument log x = 0.000(0.001)2.000(0.01)3.40(0.1)5.0; 5D]. 
Leonelli gave a similar table, log x = [0.00000(0.00001)0.00104; 14D}. 

While Gauss simplified the method of presenting Addition and Subtraction Logarithms, 
there was, of course, no adequate ground for later German writers to refer to these as 
Gauss’s, rather than Leonelli’s Logarithms. Gauss suggested that his five-place table be 
enlarged to a seven-place table, “10 to 100 times more extensive.” Such a table was published 
by E. A. MaTTHIEssEN? in 1817 or 1818, but the arrangement of the table was very bad and 
it contained numerous errors. 

When Peter Gray was preparing his six-place 
Gauss’s Tables, re-computed and re-arranged . . . [Private and temporary issue], London, 
1848, published in the following year in his Tables and Formulae for the Computation of Life 
Contingencies, his original calculations had thoroughly checked the inadequacies of Mat- 
thiessen. This work of Gray was materially extended in the second edition of Tables and 
Formulae . . ., London, 1870. He was the first to perceive the utility of Leonelli’s logarithms 

“in the calculation of life contingencies, and to him is due their introduction as well as the 
calculation of the necessary tables, which it is evident are valuable mathematically, cguet 
from the particular subject for which they were undertaken.” 

Comrie indicates the methods for applying a table of Addition and Subtraction Log- 
arithms to finding latitude and longitude, (a? — 5*)'/*, (a? + 6%)", and the solution of 
Kepler’s equation. It may also be applied to many problems in trigonometry, surveying, 
and engineering. 
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And finally, among works in English, a reference may be given to one of the numerous 
editions of 
C. BremikEr. Tables of the Common Logarithms of Numbers . . ., London, Nutt, 1887, with 
“Gauss’s Logarithms for Addition and Subtraction,” p. 535-593. Of the three types illus- 
trated by tables of Leonelli-Causs, Bremiker or Wittstein, and Zech, that of Zech is the best. 

For at least fifteen years J. T. PETERs has had a completed manuscript of Eight-place 
Addition and Subtraction Logarithms, but the opportunity for publication has not been 
offered. 

Information about a number of other tables of the types here considered may be found 
in J. W. L. Glaisher’s Report . . . on Mathematical Tables, 1873. 

a 


1 There was a second edition of this work, Supplément Logarithmique par Leonelli; 
précédé d’une notice sur l’auteur par J. HOWEL, Paris 1875, x, 75 p. This was reprinted from 
Publications de la Société des Bibliophiles de Guyenne. It was translated into German: 

o Leonelli’s logarithmische Supplemente . . . aus dem Franzésischen nebst einigen Zusdtzen 
von G. W. LEONHARDI . . ., Dresden, 1806, 88 p. There was also an Italian edition, 
Teorica e Pratica dei Logaritmi di Addizione e di Sottrazione . . . esposta con Tavole Calcolate 
a sette decimal: dall’ ingegnere, PlETRO CAMINATI, Novara, 1879, 79 p. This pamphlet was 
reprinted from o Rivista di Matematica Elementare, s. 2, ¥. 1, 1879. 

2C. F. Gauss, “‘Tafel zur bequemern Berechnung des Logarithmen der Summe oder 
Differenz zweier Gréssen, welche selbst nur durch ihre Logarithmen gegeben sind,’’ Monat- 
liche Correspondenz (Zach), v. 26, 1812, p. 498-528. The table, which fills 27 small octavo 
pages, was constructed “many years ago,” suggesting 1802 rather than 1806 when the 
German edition was published. The text, but not the table, was reprinted in Gauss, Werke, 
v. 3, Géttingen, 1866 and 1876, p. 244-246. Gauss’s table is exactly reprinted in J. Pas- 
quich, Abgekiirtzte logarithmisch-trigonometrische Tafeln ..., Leipzig, 1817; see Gauss, 
Werke, v. 3, p. 246-250. It was also added to J. DE La Lanpe’s Tables de Logarithmes a 
in an edition edited by H. G. Kohler in 0 1827; t ird ed., Leipzig, 1849, p. 213-250. 8.8 x12. 8 
cm. The text (p. 251-260) exhibits applications of such a table. 

3 Tafel zur bequemern Berechnung des Logarithmen der Summe oder Differenz zweyer 
Gréssen welche selbst nur durch thre Logurithmen gegeben sind, Altona 1817. The Latin title 
on the opposite page is dated 1818, while the Introduction has the date March 1815. The 
work was reviewed by Gauss; see his Werke, v. 3, p. 250-252. 


117[N].—FINANCIAL PUBLISHING CoMPANY, Monthly Bond Values showing 
accurate monthly values on bonds and other redeemable securities paying 
interest semiannually. For Quarter Coupon Rates from 1% to 6%. Second 
ed., Boston, Financial Publ. Co.; London, George Routledge & Sons, 
Ltd., 1941, xx, 1156 p. 25 X 36.7 cm. $150.00. 


It is fifteen years since this great work was originally published. To the first edition a 
supplement was added in 1937. The decline in levels of coupon rates and yield rates made 
necessary the appreciably enlarged second edition. 

The table now occupying pages 1-925 is one from which we can at once read off the 
monthly values of bonds paying interest semiannually, with varying Coupon Rates, Yields, 
and Periods of Maturity. The Coupon Rates = 1(i/4)6%. The Yields for coupons (a) 
1 — 22% = .00(.05)5.70%; (b) 3 —44% = 1.00(.05)6%, also 8ths 2.50 -—6%; (c) 
44 — 6% = 1.50(.05)63%, also 8ths 3 — 64%. The Maturities vary from 1 month to 40 
years and 3 months by months, from 404 years to 50 years semiannually, from 50 years to 
65 years annually, and from 70 years to 110 years by 5 years. As it is unusual to find many 
bonds of more than 40 year maturity this covers most practical needs. Ample explanation 
is given in the introductory pages to cover any maturities beyond 40 years up to the 110 
year limit. The table is to six places of decimals throughout. Convenient thumb indexes 
enable one to turn immediately to the pages with the desired coupon rate. 

On every page there are 12 columns (one for each month) and 115 lines making 1380 
entries. But these are so spaced in groups of fives that the page is singularly easy to use in 
practice. The columns of entries are facsimilies of type-script, to which many checks had been 
applied, and errors in type-setting being avoided in this way, the number of errors already 
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reported is exceedingly small; see MTE 18. In black-face type, on each side of every page 
and every line, are the 115 yields considered. At the bottom of each page are given under 
appropriate columns the obviously desirable amounts of accrued interest on $100 for 
1, 2, ---, 5 months. 

Pages 926-1108 are really part of the preceding table, filling in the gaps in the bond 
values in the lower yields [= .00(.05)1.00 or 1.50] for coupon rates from 3% to 6%, at the 
same intervals of time. 

Following this (p. 1109-1153) is a set of tables for non-interest bearing securities (0%) 
to show their values at yields = 0(.05)5.70%, intervals for the same time period as in the 
preceding tables. And finally there is a ‘‘Table of factors” (p. 1154-1156), to 8D, to be used 
in computing values for one to thirty days between the even months for yields from 0% 
to 5.7% at .05% intervals. 

The manner of preparing the tables was almost wholly dictated by trade practice and by 
what would be most intelligible to bond dealers. The method of computation involved the ~ 
addition of the constant difference between values at successive coupon rates. This constant 
difference was computed to a minimum of two decimal places beyond the six decimal places 
contained in the table, and in most cases to more than two extra decimals. The mathe- 
matician notes with interest that trade practice calls for the addition to every invoice, of 
simple interest at the coupon rate. Hence it became necessary to compute the theoretically 
correct ‘“‘flat’’ value and from this to subtract simple accrued interest. Thus when an invoice 
is made for bonds sold at a yield basis the sum which changes hands will be the theoretically 
correct amount, but it will be itemized in two distinct elements neither of which may be 
considered mathematically exact. This feature is most readily observed in the case of any 
bond selling to yield its coupon rate. The value is shown as 100 on coupon dates only, and 
for intermediate dates at something less than 100; but if the accrued interest at the bottom 
of the column be added a correct result will be obtained. 

Another curious feature is that in spite of solicitude for correct semiannual compound 
interest in the body of the book, values for less than six months’ maturity are shown at true 
discount, in accord with trade custom of dealers and banks. Again, the editors believed that 
labelling a present-worth table as the value of a 0% bond evokes a clearer concept to the 
average bond dealer, than a proper labelling would. 

Without a volume of this character, dealing in municipal securities would be a somewhat 
difficult process. For example, a dealer offers a 33% coupon municipal having 20 years and 
4 months to run on a 1.95% yield basis. What dollar price must the purchaser pay? 
Turning to the 33% coupon table, find the page “20 years, column 4 months and on the 
margin the yield 1.95%” and the price is given as 125.913848 or $1259.14 per $1000 bond. 
Also at the bottom of the column is given the accrued interest in monthly intervals, $.583333 
per $100 or $5.83 per $1000 bond. The methods of adjusting for fractional parts of a month, 
and for many other special cases of practice, are explained clearly in the introductory pages. 
These pages contain also quite explicit instructions involving the use of adding or accounting 
machines where these are available. 

This volume of tables, which is the most extensive of the kind, is of value particularly 
to dealers and those handling state and municipal securities where price is usually quoted 
ona yield basis, also for insurance companies where exact computations of yield are essential. 
The tables are obviously of importance also to the actuary or accountant in making up 
amortization schedules. The volume is strongly bound in full pig skin, and lies open easily 
to the table desired. It is too bulky for use outside of an office or a library. For the con- 
venience of salesmen and others there is Pocket Edition of Monthly Bond Values with tables 
of coupon values from 4% to 5% at }% intervals. These tables are in condensed form and 
show values to two decimal places only and cover values at monthly intervals up to 12 
years for coupon rates of 4% and #% and for the remaining coupon values at monthly 
intervals up to 28 years, quarterly 28-30 years, and semiannually 30-59 years. 

F. J. Hunt 
108 Anthony St., 
Providence, R. I. 
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118[L].—Zax1 Murs, Tables of Legendre Associated Functions, (Fouad I 
University, Faculty of Science, no. 4). [Cairo], E. & R. Schindler, 1941, 
xii, 286 p. 22.8 X 29.5 cm. Price unknown. The Egyptian University at 
the Zaafaren Palace in Cairo, was founded as a free university in 1908 
and adopted as a state university in 1925. 


Out of Egypt in the midst of war and issued at a time when the enemy was not far from 
the gates of Cairo, there comes a set of unusual tables, values of the Legendrian associated 
functions, P,™(x), n, m = 1, 2, 3, ---, 10, over the range x = 0.000(0.001)1.000. Before 
describing in detail the contents of this volume it might be well to recall some of the historical 
facts relating to these functions. 

The polynomials, P,(x), associated with the name of A. M. LEGENDRE (1752-1833), 
were first employed in 1784 by LEGENDRE, who applied them to the determination of the 
attraction of solids of revolution, in ‘Recherches sur l’attraction des sphéroides homogénes,”’ 
Mémoires des Savants Etrangers, v. 10, 1785. Inspired by this memoir LAPLACE (1749-1827) 
published in this same year one of the most remarkable memoirs ever written. It is here that 
Legendre’s spherical harmonics enter into his discussion and Laplace’s equation first ap- 
pears, though only in polar coordinates. Legendre’s name is thus properly ascribed to the 
new harmonic functions; nevertheless some writers still refer to the polynomials as Laplace’s 
coefficients. 

The functions, P,,(x), are the polynomial solutions of Legendre’s differential equation 


dx? 
which arises in connection with the solution of Laplace’s equation in spherical coordinates. 
For this reason the solutions of the equation are frequently referred to as surface zonal 
harmonics, or spherical functions. The second solution of Legendre’s equation is usually 
designated by the symbol Q,(x), and is called a zonal harmonic of second kind. 

The Legendrian associated functions of mth order and mth degree are the functions 

'x™(x), which are derived from the Legendrian polynomials by the relationship: 


(1 — x*) + 0, 


Associated functions of second kind are designated by the symbol Q,(x), and are 
defined by the formula just written down where P,(x) is replaced by Q,(x). 

The associated functions of both kinds are sometimes designated by the symbols 
Pa,m(x) and Qn,m(x), although the symbols used above are the more common. 

In some applications it is more convenient to replace x by cos @, and extensive tables 
have been computed for the Legendrian functions in their trigonometric form. When @ is 
the variable, the functions are customarily represented by the symbols P,(cos 6), Qn(cos 6), 
P,™(cos 6), and Q,"(cos 6). 

In addition to their use in problems of attraction, the Legendrian polynomials have had 
an extensive application in numerical quadrature, a use first indicated by K. F. Gauss. In 
this application the zeros of P,(x) become of great significance and considerable attention 
has been devoted to their computation. Compare RMT 92. 

The volume under review provides tables of the associated functions, P,™(x), for values 
of m and m = 10 over the range x = 0.000(0.001)1.000 to an accuracy which varies from 8 
decimal places in the beginning of the table to 3 decimal places at the end of the table. 
For those who may contemplate using the tabulated values, the following schedule shows 
the number of decimal places, D, corresponding to values of m and n: 


m =1,n = 1(1)10,D = 8.m =2,n = 2(1)9, D = 8;n = 10,D = 7.m = 3,n = 3(1)6, 
_D=8; =7(1)9, D=7; = 10, D=6. m=4, n = 4(1)6, D = 8; n = 7, 8, D =7; 
n=9,D=6;n=10,D=5.m=5,n=5,D =8;n =6,7,D=7;n =8, D=6; 
n=9,10,D =4.m=6,n =6,D =8;n =7,D =7;n = 8, D=6;n =9,10,D =4. 
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m=7,n=7,D =7;n =8, D=6;n =9, 10, D=4.m=8,n =8 D=6;2 =9, 
D=4;n =10,D=3.m=9,n =9, D =4.m = 10,2 = 10, D = 3. 
The usefulness of the table is considerably enhanced by the inclusion of values of the 
modified second difference 
M” = A” — 0.1844, 


where A” and A® are the second and fourth central differences. The theory of this useful 
difference has been set forth by L./J. Comrie in “Interpolation and allied tables” (reprinted 
from Nautical Almanac for 1937), London, H. M. Stationery Office, 1936, p. 928. 
The modified differences are used in Everett’s interpolation formula as follows: 
S(%o + €/1000) = (1 — E)f(xo) + Ef(xs) — Mi” — w(t — Me”, 
where x» and x; are consecutive values of x, and where u(t) = 2£(1 — #2). 

The tables contain (pp. 282-283) values of u(t) to 7 decimal places over the range™ 
= 0.000(0.001) 1.000. There is no indication as to the source of these values, although 
and other functions used in Everett’s formula when higher differences are employed, have 
been tabulated by A. J. Thompson in Tracts for Computers, no. V, 1921; second ed., 1943. 

In the construction of the tables, the values corresponding tom = m = landn = m = 2 
were first computed from the formula 


2m)! 
Pa™(x) = (1 — x4). 


Then, for other equal values, there was used 
= (2m + 1)(2m + 3)(1 — 2*)Pu™(x). 
In the computation of the functions for the other cases where m > m, there was em- 
ployed the following iteration formula: 
2n 
n—m 


which, when n = m + 1, reduces to 
= (2m + 1)xPa™(x). 


With respect to the computations the author makes the following comments: 

“In view of the fact that the calculation depends on the repeated use of the recurrence 
formulae involving the accumulation of small errors, it has been found sufficient to calculate 
P(x) = (1 — x*)# to 16 decimals and P;*(x) exactly in order to obtain the other functions 
accurate to 8 decimals except those for which » = 9 or 10, m = 5. The latter are accurate 
to 4 decimals only, the values of the entries being large.” 

The author believes that the accuracy of his tables is high. He states that the original 
calculations were submitted to L. J. Comriz, who differenced 6,000 entries. Of the ten errors 
thus revealed, nine were due to transcription. The values were also compared with those 
computed by H. Tauiovist over the range 0.00(0.01)1.00. This comparison ‘shows agree- 
ment as far as his figures go,”” but see MTE 17. All of the calculations were carried through 
twice either by the author, or by others. 

In view of recent interest in Legendrian functions it may not be out of place to append 
a short account of the origin and present status of tables relating to them. See also UMT 14. 

The first extensive set of values of P,(x) was computed by J. W. L. GLaIsHER and 
published in the Report of the British Association for the Advancement of Science in 
1879. This table gives the complete values for P,(x), m = 1, 2, ---, 7, over the range 
x = 0.00(0.01)1.00. Four-decimal abbreviations are to be found in JAHNKE and EmpE, 
Funktionentafeln, first ed., 1909, p. 83-84; second ed., 1933, p. 185-186; third to fifth ed., 
p. 119-120, as well as in W. E. Byerty, An Elementary Treatise on Fourier’s Series . . ., 
1895, p. 280-281. Five-decimal abbreviations are given by J. B. DALE in Five-Figure Tables 
of Mathematical Functions . . ., London, 1903, p. 78-79. 
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The most extensive calculations of the Legendrian functions are due to H. TALLovist, 
who, in addition to a treatise: Grunderna af Teorin for Sferiska Funktioner, jimte Anvand- 
ningar inom Fysiken, Helsingfors, 1905, xv, 436 p., which is rich in special tables, has 
produced the following monographs: 

(1) “Tafeln der Kugelfunctionen P,(x) und ihrer abgeleiteten Functionen,” Acta 
Societatis Scientiarum Fennicae, v. 32, no. 6, 1906, 27 p. 

(2) “Tafeln der Kugelfunktionen P,(cos @),” ibid., v. 33, no. 4, 1908. 

(3) “Tafeln der abgeleiteten und zugeordneten Kugelfunktionen erster Art,’ ibid., 
v. 33, no. 9, 1908. 

(4) “Tafel der 24 ersten Kugelfunktionen P,(cos @),” Societas Scientiarum Fennica, 
Commentationes Physico-Mathematicae, v. 6, no. 3, 1932, 11 p. 

(5) “Tafel der Kugelfunktionen P25(cos @) bis P32(cos @),” ibid., v. 6, no. 10, 1932, 5 p. 

(6) “Sechsstellige Tafeln der 16 ersten Kugelfunktionen P,(x),” Acta Societatis 
Scientiarum Fennicae, n.s. A, v. 2, no. 4, 1937, 43 p. 

(7) “Sechsstellige Tafeln der 32 ersten Kugelfunktionen P,(cos @),” Acta Societatis 
Scientiarum Fennicae, n.s. A, v. 2, no. 11, 1938, 43 p. 

In monograph (1) TaLLovist added to the computations of GLAISHER the exact values 
of P;(x) over the range x = 0.00(0.01)1.00. He also tabulated the derivatives of P,(x), 
namely P,,™(x), for m and n from 1 to 8 exactly over the same range. In monograph (6) he 
gave six-place values of P,(x) for m = 1(1)16 over the range x = 0.000(0.001) 1.000. 

The values of P,(cos @) were computed by A. LopGe in “Values of Legendre’s func- 
tions from Po to P2 at intervals of 5 degrees,” R. So. London, Phil. Trans., v. 203A, 
1904, p. 100-101, to 7 decimal places. These were later extended to 10 decimal places, 
nm = 1(1)8, over the range @ = 0°(1°)90° by TaLLovist (2). To these were added the compu- 
tations to m = 32, items (4) and (5), and in 1937 there appeared the six-place tables (7) 
computed at intervals of 10’. We also note a table for m = 1(1)7 computed to 4 decimal places 
at intervals of 1° by C. E. HoLianp, P. R. Jongs, and C. G. Lams, Phil. Mag., (5), v. 32, 
1891, p. 512-523; also, Phys. So., London, Proc., v. 11, 1892, p. 221-233. This table, accord- 
ing to TALLQVIST, contains “many and not wholly unimportant errors.’’ It was reprinted 
by BYERLY, p. 278-279, and was the first table of this type. 

To TALLoQvist we are also indebted for the derivatives P,™(cos 0) for m and n from 1 
to 8 to 10 and 11 significant figures over the range @ = 0°(1°)90°. The computation of the 
first derivative was first given by C. C. Farr, R. So. London, Proc., v. 64, 1898-99, p. 192- 
202, to 4 decimal places for m = 1 to 7 and to 3 decimal places for @ = 1° to 90°. This table is 
reprinted by JAHNKE and EmpE, first ed., p. 88-89; second ed., p. 190-191. Since we have 
4P,,(cos 0)/d@ = — P,1(cos 6), this table is also included in tables of the associated functions. 

The values of the associated functions, P,”(x), were first given by TALLQvist in his 
treatise mentioned above, and later in monograph (3) over the range x = 0.00(0.01)1.00. 
‘For even orders (m) these values are either exact or computed to 10 or 11 significant figures; 
for odd orders they are given to 10 or 11 significant figures. A companion table for P,(cos 6) 
gives values to 10 or 11 significant figures at intervals of 1° for values of m and n from 1 to 8. 

It is unfortunate that Tallqvist’s treatise mentioned above is not more readily available 
in American libraries, since it contains 50 pages of tables, some of which are not found 
elsewhere. In addition to abbreviated tables of the polynomials and the associated functions 
in both forms to 4 decimal places through m = 8, and the derivatives of P,(x) through 
n = 8, there appear tables of the roots of the derivatives and of P,(x) to 7 decimal places 
through n = 8. Tables are also provided to 4 decimal places over the range x = 0.00(0.01)1.00 
of the conjugate zonal harmonic of second kind, Q,(x), and its companion function Z,(x), 
nm = 1(1)8, defined by the equations 


Que) = Pa(2) log 


-5 


Ste =) 


| 
| 
i 


whee 


RECENT MATHEMATICAL TABLES 119 


In a work entitled Tables de Fonctions Sphériques et de leurs Intégrales, Paris, 1933, xxx, 
95, 160* p., G. Pr&vost provided tables of integrals of the Legendrian polynomials and 
the associated functions, as well as values of the functions themselves, over the range 
x = 0.00(0.01)1.00, to from 5 to 6 significant figures, for m and m = 1(1)8. 

Very extensive calculations of the zeros of P,(x) have been made, a partial bibliography 
of these being given in this MTAC, p. 51-52. The most conspicuous omission is the important 
computations of E. R. Smit, “Zeros of the Legendre polynomials,” Iowa State College, 
Jn. Sci., v. 12, 1938, p. 263-274, where 6-place values are given through m = 40. Interesting 
new asymptotic approximations were used in the determination of these values. 

The most extensive collection of values pertaining to the Legendrian functions is 
contained in an unpublished manuscript in excess of 400 pages in the possession of the writer 
of this review, which, in addition to the tables mentioned above, contains the following 
new computations: 7 

Values of P,(x) and P,(cos 6), m = 1(1)8, to 10D over the range x = 0.000(0.001) 1.000, 
together with first and second differences. Values of the integrals I,(x) = Jo*P,(t)dt, and 

n(x) = Se P,™(t)dt, m, n = 1(1)8, mainly to 10D, but in every case in excess of 10 sig- 
nificant figures, over the range x = 0.00(0.01)1.00, the derivatives P,™(cos @), to 6 or 
more significant figures, over the range @ = 0°0’(10’)90°, and P,™(cos @), to 6 or more 
significant figures over the same range, m, m = 1(1)6. The last two tables are from an 
unpublished manuscript of TALLovist kindly furnished the writer of this review. 

Since. the heroic computations of TALLQvist have secured for him a permanent place 
in the history of zonal harmonics, it might not be out of place to include a short notice about 
his work. AxEL HENRIK HJALMAR TALLQVvIsT was born February 21, 1870. He entered the 
University of Helsingfors in 1885, and became Magister and Doctor Phil. in 1890. He was 
appointed Professor of Mechanics at the Polytechnicum in 1891, Professor of Physics at 
the University of Helsingfors, 1907, and Promotor, in 1935. His first publication was a 
two-volume text on technical mechanics, published in 1895 and 1896. His bibliography 
contains more than 200 items covering many branches of mathematics and mathematical 
physics. His treatise on spherical harmonics, referred to above, is an outstanding work on 
the subject with many original features. Since the appearance of this book in 1905, TALL- 
Qvist has labored long, more than a third of a century, to supply the tables needed in the 
application of zonal harmonics to mathematical physics. His contribution is certainly one 
of the most extensive calculations ever undertaken by a single individual. 


119[E, L, S].—J. SHERMAN and R. B. Ewe t, “A six-place table of the Ein- 
stein functions,” Jn. Phys. Chemistry, v. 46, 1942, p. 641-662. 


It is stated that the calculation of thermodynamic quantities from spectroscopic data 
has become a recognized method for determining the energy, heat capacity, entropy, and 
similar properties of simple molecules. In this connection Einstein functions are extensively 
used in calculating vibrational contributions. These functions are E = x/(e* — 1), e*E’, 
and — In (1 — e*) and they are tabulated for x(= huo/kT) = (0(0.005)3(0.01)8(0.05)15; 
6D]. Brief tables were available earlier in Landolt-Bérnstein, Physikalisch-Chemische 
Tabellen, 5th ed., 1st Suppl., Berlin, 1927, p. 702-707; and W. Nernst, The New Heat 
Theorem: its Foundations in Theory and Experiment, transl. G. Barr, New York, 1926, p. 
246-251. In the construction of the table use was made of F. W. Newman's tables of the 
exponential functions (1883, 1889) and of three different expansions for the logarithm. A 
summary of the application of the table to the calculation of thermodynamic quantities is 
given by L. S. Kassel, in his ‘The calculation of thermodynamic functions from spectroscopic 
data,” Chemical Reviews, v, 18, 1936, p. 277-313. 


H. B. 
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120[L, S].—M. Korant, A. AMemrya and T. I. ‘‘Tables of integrals 
useful for the calculations of molecular energies,”” Phys. Math. So., Japan, 
Proc., s. 3, v. 20, 1938, extra no. 1, 70 p. II. idem., v. 22, 1940, extra no. 1, 
28 p. 18.8 X 25.8 cm. 


To some extent the Tables of I form a continuation of the work of S. Ikehara who com- 
piled the tables of the functions of Zener and Guillemin ! used in the work of Rosen.? The 
tables were used by T. Invi in his work * on the hydrogen molecule and he cooperated with 
the authors. He says that the integrals occurring in his work can be reduced to some basic 
integrals tabulated by Bartlett, Furry and Rosen ‘ but their tables are not sufficiently ex- 
tensive. Values for the functions A,(x), B,(x) in which exp (— xt)i* is integrated over 
the ranges, 1, © and —1, 1 respectively, are given for x = [0.25(0.25)8; 10D] this range being 
extended in the first case by 8.0(0.5)13(1)16(2)20(1)22, 24. These tables go far beyond the 
tables with 6S used by Rosen though these latter go up to » = 20 instead of m = 15 but 
for n > 10 are for x = 1(1)10, 12, 14 only. In the 4-place tables of A. Unséld § 2 has only 
the values 1, 2, 3 and x = 0(.1).2, .5(.5)2, 3. In the work of R. D. Misra * 2n = — 9(1)6 
and x ranges from x to 47, 7D being given while in a second table 2n = — 1(1)6 and 4x 
runs from 3x to 11x. Values of A,(x) form = — 3, — 2, — 1 were given by E. Gold 7 in 1909. 

The paper also contains a table of the exponential integral — Ei(— x), for x = [0.5- 
(0.5)18(1)23, 25; 18D], and other values of x, 24-48, p. 15 of II. The integrals involving 
Legendre functions, Table 4, are given, for x = [0.5(0.5)1(0.25)7; 8D]. In the tables of double 
or triple entry, 5—6, the order of the Legendre functions ranges from 0 to an integral value 
not exceeding 8 while the other integers have no more than 4 values. In the Tables of II 
13 integrals for hydrides are tabulated. P. 16-28 are occupied with lists of errata in I. 

H. B. 
1C. ZENER and V. GuILLEMIn, “The B-state of the hydrogen molecule,” Phys. Rev., 
s. 2, v. 34, 1929, p. 999-1009. 

2N. ROosEN, “Normal state of the hydrogen molecule,’’ Phys. Rev., s. 2, v. 38, 1931, 
p. 2099-2114 (see especially p. 2109-2110). 

3 T. Invi, “A contribution to the theory of the hydrogen molecule,’’ Phys. Math. So. 
Japan, Proc., s. 3, v. 20, 1938, p. 770-779. 

. H. Bartiett and W. H. Furry, “Valence forces in lithium and beryllium,” Phys. 
Rev., s. 2, v. 38, 1931, p. 1615-1622. 

5A. UNs6LD, “Zur Deutung der Intensitatsverteilung in den Fraunhoferschen Linien, II. 
Teil: Die Intensitat der Linienfliigel-approximative Lésung der Schwarzschildschen In- 
—os fiir eine beliebig geschichtete Atmosphire,”’ Zeit. f. Astrophysik, v. 4, 1932, 
” ®R. D. Misra, “On the stability of.crystal lattices,’’ Cambridge Phil. So., Proc., v. 36, 
1940, p. 173-182. 

E. Gotp, “The isothermal layer of the atmosphere and atmospheric radiation,” R. So. 
sialon, Proc., 82A, 1909, p. 62. 


MATHEMATICAL TABLES—ERRATA 


On the pages indicated Errata are listed for tables by the following: 
J. R. Atrey (p. 72), A. Gray and G. B. MaTtHews (p. 74), D. F. E. MEIssEL 
(p. 74), A. G. WEBSTER (p. 70), and R. W. WILLSON and B. O. PEIRcE (p. 74). 
In this issue we have referred to Errata in RMT 113 (JAHNKE and Empe). 


13. L. M. MILNE-THomson, Standard Table of Square Roots. The square 
roots to eight significant figures of all four-figure numbers, with printed 
differences, London, Bell, 1929. 

The following error was discovered by Mr. T. WHITWELL: 
Difference following 239.7, for 3329, read 3229. 

This error was found because Mr. Whitwell uses a precaution that I have always advocated 

when doing linear interpolation with a calculating machine. The preceding tabular value is 
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set, then the difference, which is added (or subtracted) and the number in the product 
register compared with the next tabular value. Any error of setting, or in the difference, or 
in the direction of turning, is immediately revealed by the lack of agreement between the 
machine and table. The interpolation is, of course, completed by changing the decimal part 
of the multiplier from ciphers to the required interpolating factor. 

24 May 1943 L.j..c. 


14. J. T. Peters, Sechsstellige Tafel der trigonometrischen Funktionen, ---, 
Berlin, 1929. 


In April 1939 Mr. Peters reported that a second edition of this volume was about to 
be published, and asked me to send him a list of errors which I had noted. I sent him at 
once the following: 

P.49 Sec. 4°10’30”, for 2261 read 2661 
P.55 Argument following 5°14’20”, for 50’ read 30” 
P.61 Sec 6°17'00”, for 0.006043 read 1.006043 
P. 129 Sec 17°39’50”, for 8464 read 9464 
P. 234 Headings cosec and cotg to be transposed. 
L.j.c. 


15. H. M. Parxuurst, Astronomical Tables, comprising Logarithms from 3 
to 100 decimal places . . . Revised edition, 1876; see RMT 86. 


It has been already noted (p. 20) that in this volume Parkhurst gives log N, to 102D, 
for 90 values of N = 109. Although the results given are to 102D both Tables III and IX 
(+II and XIII) in which they occur are headed “Logarithms to 100 decimal places.” Thus 
accuracy to more than 100D is not claimed. I can now report concerning the exactness of 
the printed results for everything except the digits in the 101st and 102d decimal places for 
the logarithms of 73, 89. It was found that to 100D all of the values of log N, except 3, 
were entirely correct and also that in the following 9 of these cases the results were correct 
to 102D: 7, 9, 11, 21, 33, 63, 67, 74, 94. For values of N greater than 19 the first 15 figures 
of log N are given in Parkhurst’s Table XIII. The fact that these figures which are omitted 
from Tables III and IX are all correct was established by the Editor who compared the 
data with the Peters-Sharp logs. (A reproduction of Table XIII was not immediately avail- 
able to me.) We tested the 20-place Table II independently [NV = 1, 2, ---, 19], and detected 
no error. 

The three non-terminal errors in log N were as follows: In log 51 the 9 in the 82d 
decimal place should read 8, the correct environment being ... 88 386 .... In log 83 there 
is an error of +1 in the 91st decimal place which affects three consecutive figures. Park- 
hurst’s “‘... 18140 01079 ...”" should read ... 18139 91079 .... In log 91 Parkhurst’s 3 
in the 57th decimal place should read 4 since the correct environment is ... 04825 .... 

In order to correct the values of the 79 remaining log N to 102D, the digits after the 
100th have to be decreased or increased by amount (—15 to +15) for the N’s which are 
indicated as follows: (—15) 71; (—11) 82; (—10) 41; (—9) 39; (—8) 97; (—7) 52, 76, 78; 
(—6) 26, 38, 57, 64, 91; (—5) 13, 32, 48, 59, 65, 86; (—4) 19, 43, 58, 72, 106, 109; (—3) 8, 
16, 24, 36, 56, 87, 88, 95; (—2) 4, 6, 12, 29, 44, 53, 54, 103; (—1) 2, 3, 14, 18, 22, 27, 28, 
42, 66, 68, 81, 84, 107; (+1) 5, 15, 17, 34, 37, 45, 47, 51, 62, 75, 77, 105; (+2) 25, 35, 49, 
55, 79, 83, 93; (+3) 31, 46, 85; (+4) 69; (+5) 23; (+15) 61. 

The following account of a device for determining simultaneously the number and 
decimal positions of all the non-terminal errors in a table of logarithms of integers may be 
interesting because it was new to me and it is very effective as compared with the more 
obvious but crude scheme of computing a complete table of higher approximation in order 
to test a table having a smaller number of decimal places. 

The basic idea is that 


k 
log = + j)] = log [(m + k)!] — log [n!]. 
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The device consists in adding the logarithms of the table under investigation and then 
comparing the sum with standard values of the logarithms of the factorials. At the be- 
ginning of a table = 0. In the special case of Parkhurst’s Tables III and IX the values of 
log N are omitted for N equal to all multiples of 10 and for N = 92, 96, 98, 99. The missing 
logarithms were supplied frém my own master values after conversion to the base 10. 
More specifically, the characterizing value of log (100!) was computed by the formula 
12 + 14 log 2 + 7 log 3 + 3 log 7 + log 11 + log 23 + (sum of Parkhurst’s log’s inclusive 
of log 97). Comparison of this sum with the converted value of my published 153-place 
approximation to log (100!), as obtained from Stirling’s series for log '(x), showed only 
three discrepancies. As indicated earlier the numerical differences had the common value 
1 and they occurred respectively in the 57th, 82d and 91st decimal columns. The above 
comparison also showed that the algebraic sum of the two undetermined terminal errors 
equals +-15.1 in terms of unity in the 102d decimal place. 

The determination of the second codrdinate of the separate errors was effected by com- 
parison of Parkhurst’s logarithms with juxtaposable photostatic copies of the 61-place 
Peters-Sharp table, of A. Grimpen’s 84-place table for prime N’s, and with my own multi- 
fariously checked log’s. For values of N greater than 100 the complete comparisons were 
made of necessity with my own converted log’s. 

The omnibus exploratory test explained above is a necessary but not a sufficient condi- 
tion for perfection in a table since the theoretical probability of the miraculous occurrence 
of mutually compensating non-terminal errors is not zero. This qualification applies to the 
present concrete case of Parkhurst’s tables only for the logarithms of composite values of NV 
in the range of figures from the 60th (end of Sharp’s table) to the 79th decimal place, a 
region within which it was not considered necessary to attain purely theoretical sufficiency 
by multiplying every figure of my natural logarithms by the common modulus. 

H. S. UBLER 
Yale University 


FINAL NoTE. On June 9, 1943 I finished calculating to about 121D the natural logarithms 
of 73 and 89, and in converting these data to the base 10 in an interval including the 100th 
dec. place. It was found that the respective pentads beginning at the 101st place are 10647 
and 63094, hence Parkhurst’s terminal pairs ‘‘23’” and ‘‘66” are too high by about 12 and 3 
in terms of unity in the 102nd place of decimals. Finally as an omnibus check on all of my 
base 10 sequences, which had been used as standards of comparison for Parkhurst’s ter- 
minal figures, I computed from them a section of log (100!) on the basis of the exponents of 
the 25(>1) prime factors of 100!. The section thus obtained, which begins at the 96th dec. 
place, reads . . . 16858,35310,985 . . . and it agrees absolutely with the figures in the 
approximation to log (100!) as derived earlier from my coefficients of Stirling’s series. 
H. S. UHLER 


16. PROJECT FOR THE COMPUTATION OF MATHEMATICAL TABLES, Tables of 
Probability Functions, v. 2,1942. Compare RMT 91. 
P. 46, for 0.2335, read 0.2235 
P. 63, for 0.0370, read 0.3070. 
D. H. L. 


17. A. H. H. Tatrevist, Tafeln der abgeleiteten und zugeordneten Kugel- 
functionen erster Art. Soc. Scient. Fennicae, Acta, v. 33, no. 9, 1908, 
67 p. Compare RMT 118. 


P#3(0.19), for 0.8917, read 2.8917 
P4(0.64), for 46.6000768, read 36.6000768 


‘For P;*(0.40), P;*(0.41) the respective correct values, — 118.19932808 and —121.11152590, 


are interchanged. 
Z. Murst, Tables of Legendre Associated 
Functions, Cairo, 1941, preface. 
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18. FINANCIAL PuBLIsHING Company, Monthly Bond Values . . ., second 
ed., Boston, 1941. Compare RMT 117. 


The following list of errors in this work was supplied by the Company, and is complete, 
according to their latest information, except for a few last-figure unit errors, which are of 
no consequence in a table of this kind. The Company has remarked that since computation 
is made arithmetically from one coupon to the next, a single original mistake may appear, 
as below, in several or all coupon rates. 


page coupon yield yrs. mos. for read 
3 1% 3.00 1 3 97.562848 97.563248 
5 1 AS 3 5 101.862678 101.862638 
47 14 3.00 1 3 97.866993 97.867493 
49 1} 45 3 5 102.709357 102.709307 - 
91 1 3.00 1 3 98.171138 98.171738 
93 1 AS 3 5 103.556036 103.555976 
135 12 3.00 1 3 98.475283 98.475983 
137 13 AS 3 5 104.402716 104.402646 
143 12 2.05 9 11 97.019650 97.319650 
179 2 3.00 1 3 98.779428 98.780228 
181 2 45 3 5 105.249395 105.249315 
223 23 3.00 1 3 99.083573 99.084473 
225 23 45 3 5 106.096075 106.095985 
267 2 3.00 1 3 99.387718 99.388718 
269 2 A5 3 5 106.942754 106.942654 
311 22 3.00 1 3 99.691863 99.692963 
313 2? AS 3 5 107.789433 107.789323 
442 34 2.85 0 8 100.323488 100.423488 
887 6 3.85 5 8 104.680108 114.680108 
926 3 45 3 5 108.636113 108.635993 
938 3} 45 3 5 109.482792 109.482662 
950 34 A5 3 5 110.329472 110.329332 
962 3? 45 3 5 111.176151 111.176001 
974 4 45 3 5 112.022830 112.022670 
986 4} 45 3 5 112.869510 112.869340 
999 43 A5 3 5 113.716189 113.716009 
1015 43 AS 3 5 114.562869 114.562679 
1030! 5 .90 1 7 106.420496 106.430496 
1031 5 AS 3 5 115.409548 115.409348 
1047 5 45 3 5 116.256227 116.256017 
1063 5 A5 3 § 117.102907 117.102687 
1079 52 45 3 5 117.949586 117.949356 
1095 6 . 3 5 118.796266 118.796026 
1117 0 7 10 


values at 2.00 and 2.25 trans- 
posed 


1 This particular error does not appear in all books. In the middle of the press run a 
super-efficient pressman decided to improve a blurred figure, and in doing so substituted 
as 
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19. J. P. L. Bourcuert, “Sur les intégrales Eulériennes et quelques autres 
fonctions uniformes,’’ Acta Mathematica, v. 2, 1883, p. 261-295. 


In this paper on the gamma function, among the six tables, to 16D, are two (p. 288-289) 
for coefficients in the following expansions of 1/I'(x): 


1/P(x) = x(x + 1)(1 + Bix + + But +---), and 

1/T(x) = x + Cx? + Cx? + 
These tables are reprinted in H. T. Davis, Tables of the Higher Mathematical Functions, 
Bloomington, Indiana, v. 1, 1933, p. 185. 

Besides checking these coefficients, it was thought desirable to compute the remaining 
few coefficients that are nonvanishing in the 16th decimal place. Corrections and new values 
for the B,’s are accurate to much within 3/4 of a unit in the 16th place, and for the C,’s 
(which are obtained by addition of B,’s) to much within 1} units in the 16th place. 

Errors in the last four figures of the B,’s: 


n for read difference 
7 6044 6042 —2 
11 0224 0290 +66 
12 5466 5469 +3 
13 3152 3149 -3 
18 0028 0030 +2 
19 2202 2207 +5 


Addition to B,’s: Bz; = + 13, and Ba = + 1 units in the 16th place. 
Errors in the last four figures of the C,’s: 


n for read difference 
12 7741 7807 +66 
13 4758 4821 +63 
14 2314 2320 +6 
15 8420 8417 —3 
19 3425 3427 +2 
20 7826 7823 -3 
21 6962 6968 +6 


Addition to C,’s: Cag = — 54, Cos = + 14 and Cos = + 1 units in the 16th place. 
E. Isaacson and H. E. SALzer 
Mathematical Tables Project, 
New York City 


20. R. A. Davis, Table of Natural Sines and Radians, Oakland, California, 
1941. Compare UMT 7. 


This table lists sin « in radians for x = [0(.001)1; 6D]. It was proofread against the 
Project’s manuscript and no error larger than two units in the last place was found in Davis. 
But the errors show a trend: Of the first 450 values there are 265 entries which are too small 
by a unit in the last place. The only other last place which is a unit too small occurs at .907. 
Among the last ninety entries from .910 to 1, are 34 values which are too large by a unit in 
the last place. The only other last place which is a unit too large is at .416. In addition, 
the last place was two units too small for the following arguments: .053, .055, .057, .058, 
061 — .064, .110, .116 — .119, .163, .171, .176, .178, .179, .181, .182, .185, .186, .189, .236, 
258, .353, .354, .396, and .441. Also the last place was two units too large for .998. Davis 
uses a + or — sign after the last place from .450 to .900. In that interval, we allowed hima 
unit leeway. Thus if his value was larger by a unit in the last place followed by a minus sign 
or smaller by a unit followed by a plus sign, it was not noted. Also if his last place value 
agreed with ours, we ignored the plus or minus sign that sometimes followed. Under these 


conditions no error was found from .450 to .900. 


I. I. Ruopes and H. E. Sarzer 
Mathematical Tables Project, 


New York City 
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References are made to unpublished tables of PETERs in RMT 116, of 
Davis in Q 5, and of Davis and of TALLovist in RMT 118. 


10[A].— HERBERT ELLIs SALZER and ABRAHAM HILLMAN (1918— ), The 


First 120 Factorials. Ms. in possession of the authors; a copy is in the 
Library of Brown University. 


The exact values of the first 120 factorials were obtained as follows: Direct multiplica- 
tion (checked by duplicate calculation at every step) was used to obtain the first 100 fac- 
torials. An additional check was performed by the independent calculation of every fourth 
or fifth factorial, and a manuscript was prepared which was proofread against the original. 


The value of 100! was found to agree with that of H. S. Uhler in Nat. Acad. Sci., Proc., © 


v. 28, 1942, p. 61. Then the table was extended to 120!. The largest table of factorials before 
this was given in J. T. Peters, Zehnstellige Logarithmentafeln, v. 1, Berlin, 1922, Anhang, 
p. 58, containing the first 60 factorials, all of which were found to be correct when proofread 
against this table. Finally, to close the gap of possible uncertainty, due to proofreading, the 
final manuscript was checked again by the independent computation of every value from 60! 
to 120!. Thus the authors are confident that no error exists in their manuscript. 

In L. Potin, Formules et Tables Numériques relatives aux Fonctions Circulaires, Hyper- 
boliques, Elliptiques, Paris, 1925, p. 836, are the first 50 factorials, which contain the follow- 
ing errors: 18! has an extra 0 at the end; 38! has the first group of ten numbers written as 
52302 26147 which should read 52302 26174; 45! has the first group of ten numbers written 
as 11962 22086 which should read 1 19622 22086; and 50! has an error in the fifth group of 
ten numbers written as 75415 68960 which should read 76415 68960. 

H. E. Satzer and A. HILtMAN 
Mathematical Tables Project, 
New York City 


11[D, G].—Table of Chebyshev Polynomials. Preliminary Ms. prepared by, 
and in possession of, the PROJECT FOR COMPUTATION OF MATHEMATICAL 
TABLEs, 50 Church St., New York City. 


The Chebyshev Polynomials C,(x) and S,(x) are defined by the expansions of 2 cos n@ 
and sin (m + 1)0/sin @ in terms of the argument x = 2cos#@. These polynomials were com- 
puted to twelve decimal places for m = 1, 2, ---, 12 and for x ranging from 0 to 2 at intervals 
of 0.001. The only similar table in existence is contained in an unpublished manuscript of 
the Tables Committee of the British Association for the Advancement of Science; this 
manuscript lists twelve-place values of C,(x) for m = 1, 2, ---, 12 and for x ranging from 
0 to 1 at intervals of 0.01. 

Chebyshev Polynomials arise in many phases of mathematical analysis and are of 
particular practical importance as a tool for approximating empirical or analytic functions 
and for evaluating power series. An exhaustive bibliography of Chebyshev Polynomials is 
contained in A Bibliography on Orthogonal Polynomials (Nat. Res. Council, Bull. no. 103, 
1940). See also G. SzEG6, Orthogonal Polynomials (Amer. Math. So., Colloquium Publications, 
v. 23, 1939). 


A. N. Lowan 


12[L].—Tables of Jacobi Elliptic Functions. Preliminary Ms. prepared by, 
and in possession of, the PROJECT FOR COMPUTATION OF MATHEMATICAL 
TABLEs, 50 Church Street, New York City. 


These are Tables of sn(u, k), cn(u, k) and dn(u, k), to 15D, for k* = 0(0.01)1, and for 
u = pK, where p = 0.01 and p = 0.1(0.1)1, and K is the complete elliptic integral corre- 
sponding to the modulus k. 
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From the computed values one can generate the values of the functions under considera- 
tion for imaginary arguments. Specifically, sm(iu,k) = isn(u, k’)/cn(u, k’); cn(iu, k) = 
1/cn(u, k'); dn(iu, k) = dn(u, k’)/cn(u, k’) where k’ = vi — #. It is also possible to generate 
the values of the function am(u, k) = sin sn(u, k). 

With the aid of the known addition formulae, it is further possible to generate the 
values of the function in question for all p’s ranging from 0 to 1 at intervals of 0.01. 

The three Jacobi elliptic functions were computed by two independent methods. 
(a) From their known expressions in terms of the theta functions, as defined in E. T. 
Whittaker and G. N. Watson, Modern Analysis. (b) By the AGM (arithmetico-geometric 
means) method, involving the sequence of number-triplets, as described, for instance, in 
L. V. King, On the Direct Numerical Calculation of Elliptic Functions and Integrals, Cam- 
bridge, Univ. Press, 1924. 

The following auxiliary functions were also obtained in the course of the computations. 


(1). The complete elliptic integral K for k? = [0(0.01)1; 17S]. 


(2). g® = [Exp (— +K’/K)}*" for m = 1/4, 1/2, 1, 2, 3, 4, 6, 8, 9, 12, 16, 20, and 25, to 17D 
— log g to 15D; K’ is the complete elliptic integral corresponding to the complementary 
modulus k’. 


(3). The three theta functions: 32(7u/2K, g), 33(xu/2K, g), and &(xu/2K, for u/K = 0, 
0.01, and 0.1, to 17S. 


(4). The first six sets of number triplets of the AGM methods to 17D. 


__ The only known tables of the Jacobi Elliptic functions sn, cn and dn are contained in 
L. M. Milne-Thomson, Die elliptischen Funktionen von Jacobi, Berlin, Springer, 1931. This 
is a 5-place table in which sm, cn and dn are tabulated as a function of k* ranging from 0 to 1 
at intervals of 0.1 and of the argument u whose range and interval are as follows: for the sn 
function, u = 0(0.01)3(0.1)6.5; for the cn function, u = 0(0.01)3.50; for the dn function, 
u = 0(0.01)4. 

Some of the auxiliary functions listed under (1), (2), (3), and (4) are also tabulated 
in K. Hayashi, Tafeln der Besselschen, Theta-, Kugel- und anderer Funktionen, Berlin, 
Springer, 1930. 

A. N. Lowan 


13{I, K].—Table of coefficients for inverse interpolation. Ms. prepared by, and 
in possession of the PRoyEcT FOR COMPUTATION OF MATHEMATICAL 
TABLEs. 50 Church St., New York City. 


Although many tables have been published to facilitate direct interpolation, there has 
been to date no table of coefficients to handle the more cumbersome task of inverse interpola- 
tion. To meet this need, the Project has computed the polynomial coefficients of the ratios 
of differences of various order in the formula obtained by the inversion of the Everett formula 
for direct interpolation (H. T. Davis, Tables of the Higher Mathematical Functions, v. 
1, Bloomington, Ind., 1933, p. 82-83). Specifically: The coefficients of the five fourth 
order terms were caiculated to ten decimal places at intervals of 0.001 of the argument 
m = (u — uo)/(ui — uo). Also a short table giving the values of the coefficients of the ten 
sixth order terms was computed at intervals of 0.1. It was not necessary to compute the 
coefficients of the second order terms, since these coefficients are tabulated at intervals of 
0.0001 in the table of Lagrangian interpolation coefficients to be published by the Columbia 
Press. Compare UMT 6, p. 94. The table here described will be of particular value whenever 
inverse interpolation is to be carried out in a table with a fairly large tabular interval. 

A. N. Lowan and H. E. Sauzer 
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14{L].—James NicHoLas SARMOUSAKIS (1912— _—+), Tables of associated 
Legendre functions with complex argument. Manuscript and film in pos- 
session of the author. Film copies at the Mathematical Tables Project, 
New York City, and in the Library of Brown University. 


The Tables of P,™(ix) and of Q,"(ix) are as follows: 


P,™(ix), x = 0(.1)5 for m = 1(1)5, x = 0(.1)2.5 for » = 6, x = 0(.1)1.4 for m = 7(1)9, 
x = 0.1(.1)1 for m = 10(1)12, x = 0.1(.1)0.8 for m = 13, 14, x = 0.7 for m = 15(1)18, 
and m = 0(1)n. 

Qn™(ix), x = 0(.1)14 for m = 0(1)9, x = 0.1(.1)1 for m = 10(1)12, x = 0.1(.1)0.8 for 
n = 13, 14, x = 0.7 for m = 15(1)18, and m = 0(1)(m + 1). The tables range from 4S 
to 20S, mostly between 10S and 15S. 


The computation of these tables was done only once and the sole auxiliary tables used 
in their preparation were the Tables of the First Ten Powers of the Integers from 1 to 1,000, 
New York, Project for the Computation of Mathematical Tables, 1938. A few miscellaneous 
checks have been made. As yet there has been no appropriate rounding off to a suitable 
number of significant figures. 

I have used the tables to evaluate relative dissociation constants of substituted benzoic 
acids which were computed on the basis of an electrostatic theory of substituent effects. 


J. N. SaRmMousakts 
Aberdeen Proving Grounds 


Such tables are also of use for computing potentials, flow and diffusion about disks and 
rods.—EbITor. 


MECHANICAL AIDS TO COMPUTATION 


6[Z].—Stmpson Leroy Brown (1881-  ), “A mechanical harmonic 
synthesizer-analyzer,” Franklin Institute, Jn., v. 228, 1939, p. 675-694. 


The 30-term synthesizer-analyzer described in this paper is perhaps the largest machine 
of its kind having so few harmonic components. It is 15 feet long and 7 feet high and 
weighs nearly a ton. This may be compared with the old Michelson and Stratton! 80-term 
machine which is about 3 feet long and 4 feet high. The reasons for building the machine 
on a large scale were to gain accuracy and to save labor. The driving mechanism, for ex- 
ample, is a train of 22 spur gears which, though of only commercial grade, are large enough 
so that more accuracy is obtained than with small expensive precision made gears. 

Two simple harmonic motions, out of phase by 90°, are generated at both ends of 15 
rotating shafts by Scotch crossheads. This gives 15 sine terms and 15 cosine terms in the 
corresponding Fourier series. The two fundamental motions are capable of amplitudes up 
to 16 centimeters while even the components of highest frequency may be set with ampli- 
tudes as great as 4 centimeters. Thus the machine draws curves well over a foot in overall 
width. The 30 harmonic motions are added, in the usual manner, by an endless fine chain 
passing over pulleys. 

Harmonic analysis to, say, 60 terms could be performed by pencil and paper methods, 
better still by commercial computing machines; it may also be done entirely mechanically 
with a 60-term analyzer. The author points out that a combination of the two methods is 
not only possible but indeed desirable, and that by using an accurate synthesizer with com- 
paratively few harmonic elements one can make analyses that include many harmonic 
elements by applying a little pencil and paper before and after using the analyzer. 

This may be illustrated as follows. Let y = f(x) be the subject of analysis, a periodic 
function defined empirically by 64 of its equally spaced arguments Yea). We wish 
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to find constants A, and B, so that 


y= fis) = z ( 41 cos + B, sin “= ) 


mci 


(1) res + be sin 32 <=) 


Before using the machine we prepare 60 constants 
am, bm, Gu’, bn’, (m = 1, 2, 15) 
defined in terms of the given y’s by 


= Ym — Y32-m — + 
Om = Ym + — — 
Om’ = Ym + Ys2-m + Ys24m + 
= Ym — + — Voum 


Setting the amplitudes a», b, on the machine we draw the curve represented by (1) and 
take therefrom 64 equally spaced ordinates (Yi, Yo, ---, Yes). Actually, no curve need be 
drawn and measured; the machine simply is stopped and read 64 times at equal angle 
intervals. Similarly am’ and bm’ lead to (¥Y1', V2’, ---, Yes’). The desired numbers A, and B, 
may now be found as simple combinations of the Y, Y’, and yie, ys2, Yas, You. Thus for 
j = 1,2, ---, 16 we have 


64A 93-1 = + + — 

64A 0; = + + + + 2(— + 
64B2; = — 


The author gives also a corresponding scheme for the analysis of a function defined by 120 
equally spaced ordinates. This involves the use of the machine four times. 

The paper contains 6 cuts taken from drawings made by the machine. Amplitudes as 
small as 1/10 mm. are said to be discernible. For extreme sensitivity in drawing, the tracing 
point is replaced by a beam of light sharply focused on photographic paper. As indicated 
above however, drawing curves is an activity of the machine which is of secondary impor- 
tance to the business of accurate harmonic analysis, but is of much more importance to the 
other uses to which the machine is put. These will be taken up in a later review. 

D, 8. 1. 


1A, A. MICHELSON and S. W. Stratton, “A new harmonic analyzer,” Amer. Jn. Sci., 
v. 155, or s. 4, v. 5, 1898, p. 1-13 + plate. 


7[(Z].—S. L. Brown and LisLE LorRENzO WHEELER (1904— ), “A me- 
chanical method for graphical solution of polynomials,”’ Franklin Insti- 
tute, Jn., v. 231, 1941, p. 223-243. 


This paper discusses the application of the harmonic synthesizer, described in the 
preceding review, to the problem of graphing a polynomial function of a complex variable 
with a view to obtaining the roots of the polynomial. 

The method is one suggested in 1928 by A. J. Kempner and applied by a specially 
constructed machine, designed by T. C. Fry, by the Bell Telephone Laboratories and called 
the Isograph.! 

The present machine is capable of handling equations up to the 15th degree. 

Briefly the method is as follows. 

Let 

z = r(cos 6 + isin 6) 
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be a complex variable. If 


f(x) = = 0 


is the equation under discussion, the machine is set to draw the map M in the w-plane of 
the circle 

under the transformation 


w= an (cos nO + isin n6). 
The number of times this map encircles the point w = — do is the number of roots within 
the circle |z| = 7. By choosing different values of r and inspecting the different maps, a _ 
value of r is soon found whose map passes through the point w = — do. This value of r is 
the absolute value of one root of f(x) = 0. The argument of this root is the angle @ read on 
the machine just as the tracing point reaches the point w = — do. 

To draw these maps the sine elements of the machine are used to drive the tracing point 
up and down while the cosine elements are used to drive the drawing board back and forth. 
The amplitude (for a fixed r) set on the mth crosshead is |a,|r”. If a, is not real, its argument 
is set by giving the crosshead an appropriate initial angular displacement. A dozen or more 
of these maps are illustrated in the paper. 

A second method applicable to real roots only is given. Suppose that one wishes to 
find all the real roots of f(x) = 0 that lie between —r and +r. Substitute x = r cos@ in 


k 
f(x) = 2 a,x". Since cos" @ is a linear combination of cos 8, cos (n — 2)0, ---, the poly- 
n=O 


nomial f(x) is reduced to a series of cosines which the machine can draw in one application. 
The real roots of f(x) = 0 (and indeed any other features of the graph of f(x) for — r =x =r) 
can then be read from this graph. The paper offers convincing evidence that this harmonic 
synthesizer is a powerful tool for studying polynomials and their roots. 

D. H. L. 


1 Bell Laboratories Record, v. 16, 1937, R. L. Dietzoxp, ‘The isograph—a mechanical 
root-finder,” p. 130-134; R. O. MeRcner, ‘The mechanism of the isograph,” p. 135-140. 


NOTES 


10. BrIGGs AND ViETA.—That Briggs in early years of the seventeenth 
century should have computed a “canon of sines,” and of other functions, 
for every hundredth of a degree in the quadrant, is a notable fact, and that 
he was probably mainly led to this choice of argument by a passage in 
Vieta’s Relatio Calendarii vere Gregoriani, Paris, 1600 is of special interest. 
The passage on p. 1 of his Trigonometria Britannica (1633) where Briggs 
tells us this is, in translation from the Latin, as follows: ‘“Therefore we divide 
any circumference into 360 parts which we call degrees; and each degree we 
divide sexagesimally into minutes and seconds, etc. But I, induced by the 
authority of Vieta ‘pag. 29. Calendarij Gregoriani,’ and on the advice of 
others, divide degrees by the decimal system into 100 primary parts, and 
each of these into 10 parts, of which each one is divided in the same way. 
And these parts give a much easier, and not less certain method of calcula- 
tion.” With too hasty a glance at this passage I made an incorrect state- 
ment (which might have been corrected with the change of a single word) 
on page 33, line 12; see Corrigenda, p. 100. Since Vieta’s Relatio of 1600 is 
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not to be found in any of the larger libraries of America, I turned to the 
reprint on p. 447-503 of Vieta’s Opera Mathematica, Leyden, 1646. But I 
failed to find anywhere a passage such as seemed to be referred to by Briggs. 
That this same difficulty was experienced by such a distinguished table- 
maker and authority as J. W. L. Glaisher, who had the original work of 
Vieta before him, is shown by the following note, in his handwriting, which 
I happen to possess (in several places the writing is difficult to read): 

“80 pages 400 [should be 40] double pp. as the rectos are numbered 1 to 40.” 

“I have looked through Francisci Vietae Fontenaeensis . . . Relatio 
Kalendarii vere Gregoriani . . . 1600 (Colophon-excudebat Parisiis . . .) 
without finding either on p. 29 or anywhere else anything about the decimal 
division of the degree. Without wanting to say that the book does not con- 
tain anything of the kind, it is not unlikely that the reference is a wrong one. 
Vieta always in the [this?] book seems to use when he can dec[imal] p[arts] 
viz. fra{ctions] with 100 . . . as the den[ominato]r.” 

On laying this note before Alexander Pogo, an outstanding authority 
not only on the Gregorian but also on many ancient calendars, he cleared 
up the difficulty completely in the following communication. 

R. C. A. 


Vieta does not always use fractions with 1000,000,000 as denominator. 
He uses them on p. 487 of the Leyden 1646 edition for the first time, with 
an explanation indicating that not all his readers are ready to accept them; 
he uses the same denominator, 1000,000,000, again three times on p. 489, 
once on p. 497; he uses 1000 as denominator three times on p. 494; 10000 on 
p. 497; 100,000,000,000 on p. 497. That seems to be all. 

And here is mathematical proof that leaf (folio, not page) 29 of the 
Paris 1600 edition is actually meant. There are 80 pages in the Paris 
edition, hence leaf 29 would correspond to p. 57 or p. 58 of the Paris edition. 
The Leyden edition contains the Relatio on pages 447 to 503, or on 57 pages, 
including a title page (p. 447); dividing 57 pages of the Leyden edition in 
the ratio of 57 : 80 or of 58 : 80, we find page 41 of the reprint, i.e., the 
page 446 + 41 = 487 of the Leyden edition. This is the page on which Vieta 
introduces his 20 maw days, after having given the clumsy execontades 
and 29535. As far as I am concerned, this is Briggs’s reference—fol. 29 
(recto or verso?) of the Paris 1600 edition; page 487 of the Leyden 1646 
edition can be quoted as the practically certain equivalent of fol. 29 of the 
Paris 1600 edition. 

The idea of a decimal fraction of a degree is obviously contained in 
Vieta who gives decimal fractions of a day, instead of hours, minutes, 
seconds, etc. Briggs did not say that Vieta subdivides a degree decimally. 
Briggs must have seen the Paris equivalent of p. 487 of the Leyden edition, 
where Vieta says: 

“At qualium dierum 297355 seu aliter 2952-92 est mensis unus, talium 
menses duodecim sunt dierum 354” 

A better case for the simplicity of decimal subdivisions of a day or of a 


' degree, avoiding the sexagesimal fractions, would be hard to find. You 


multiply 29.530,592,348 by 12, and you get the result. 


A. Poco 
arvard University Library 
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QUERIES 


4. Loc Loc Tasies.—In calculations of some physical phenomena 
the logarithms of logarithms are employed. Are there any tables to facilitate 
such calculations? 

H. K. HAMMER 
17 Howland Ave., Beacon, N. Y. 


EpitoriAL Note.—In E. CHAPPELL, Five-Figure Mathematical Tables, London, 1915, there 
is a log log table p. 171-207, and an anti log log table p. 209-219. There is also a log log table 
of numbers, with varying intervals, from 0.0010 to 1000, in S. P. GLazEnaP, Matematicheskie 
4 Astronomicheskie Tablitsy, Leningrad, Acad. Sci., 1932, p. 49-69; but it is obvious that 
this table is a reproduction of the corresponding table of Chappell, except for the omission 
of his tables of proportional parts. Chappell’s term “‘Lologs” heads each page. In P. F. 
oe Traité taire des Fonctions Elliptiques, Brussels, 1841, there is a table of 
log log g™, to 12-14D, p. 305-309. There are other tables of this type, as well as similar 
tables, with different argument, for example, by H. NaGaoxa and S. Saxurat, in Inst. 
Phys. Chem. Research, Tokyo, Sci. Papers, v. 2, 1922, p. 7-16. Where are other log log tables? 


5. TABLES OF N*/?,—We have had a request for three-halves powers of 
numbers to more than three places. The three-place table, on p. 23 of M. 
MERRIMAN, Mathematical Tables for Classroom Use, New York, 1915, is in 
our Library. We should appreciate receiving references to any table such as 
has been requested. 

U. S. NATIONAL BuREAU OF STANDARDS LIBRARY 

EprroriaL Note.—The three-place table referred to above is given also in American Civil 
Engineer’s Handbook, ed. by T. MERRIMAN and T. H. WiacI1n, fifth ed., New York, 1930, 

. 25. This is a table for N = [0.00(0.01)4.49; 3D]. It is an abridgement of tables in R. 

. Horton, Weir Experiments, Coefficients, and Formulas, U. S. Geological Survey, 
Water-Supply and Irrigation Paper, no. 200, 1907, p. 171-176, N = [0.000(0.001)1.490; 
mostly 5D], and N = [0.00(0.01)12; 4D]. Table 41 in U. S. Dept. of Interior, Bureau of 
Reclamation, Hydraulic and Excavation Tables, eighth ed., 1940, p. 126-132, gives N*”, for 
N = [0.000(0.001)1.499; 4D] and N = [1.50(0.01)19.99; 4S]. By making use of Barlow's 
Tables of Squares, Cubes, Square Roots, . . . of all Integer Numbers up to 12,500, edited by 
L. J. Comrie, fourth ed., London, 1941, 23% is found to be 12167; the square root is given 
as 110.304125. Hence this volume readily yields the three-halves power of integers, and 
other numbers, up to 23, to more than 3D. H. T. Davis, of Northwestern University, has 
manuscript tables of N+’ for N = [1(1)100, and 100(10)1000; 10D], with first differences. 
What other tables of N*/? are there, to more than 3 


6. AN ENGEL TaBLeE.—In 1928 the State Statistical Bureau of the 
Czechoslovakian Republic published Handweiser zur harmonischen Analyse, 
by LEo WENZEL Po tak, then professor in the German University of 
Prague, but now residing in Ireland. This appeared as part II of Prager 
Geophysikalische Studien, or of Cechoslovakische Statistik, v. 54, or s. 12, part 
10. 23.3 X 31.6 cm. It is a work of 72*, 98 p. most of it filled with important 
mathematical tables. On p. 28* is a reference to Zehnstellige Tafeln der 
Sinus, Cosinus und Tangenten fiir die dezimale Teilung des Nonagesimalgrades, 
calculated by ‘‘Hofrat E. Engel.” This was published in 1920 by the “General 
Direktion des Grundsteuerkatasters (Osterr.-Triangulierungs- und Kalkul- 
bureau).’’ My attention was first directed to this table by J. C. P. Miller 
in a footnote of his review in Math. Gazette, 1942, p. 228. In what public or 
private library may this Engel table be consuited? It is not in the larger 
libraries of America, or in the library of the U. S. Geological Survey. 

R. C. A. 
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QUERIES—REPLIES 


6. SCARCE MATHEMATICAL TABLES (Q 2; QR 5).—Mr. R. D. Brown Jr. 
of Blue Bell, Montgomery Co., Penna., has a copy of A, which he kindly 
loaned to us, so that a film might be made for the library of Brown Uni- 
versity, and checks with other tables might be effected. The title-page 
(corrected) is as follows: Calcul des Rouages| par| Approximation. | Nouvelle 
Méthode|par| ACHILLE Brocot|Horloger.| Paris| Achille Brocot, 6, rue du 
Parc-Royal | Londres| Lemaitre et Bergmann, 65, Cannon Street (West) | near 
St. Paul’s Cathedral| 1862. 97 p. 16.7 X 26.1 cm. On p. 47 is the half title: 
“Table de conversion en Décimales des Fractions ordinaires. A l’usage du 
Calcul des Rouages par Approximation. Méthode Nouvelle.” This is exactly 
the title of B which contained 51 pages, precisely the number of the pages 
47-97. This suggests that these pages of A were reprinted to form B. 

In MTAC, p. 22, it was noted that the Hiitte table gives the decimal 
equivalent, to 11D, of all proper fractions whose positive numerators and 
denominators are not greater than 100; and it was remarked ‘There seems 
to be no doubt that this is the form of Brocot’s original table.’”” We now find 
that the tables are practically identical in form except that the Brocot table 
is to 10D, instead of to 11D. The table occupies p. 47-89, and pages 91-97 
are in English, the half title being ‘‘Use of the table for the calculation of 
clock-trains by approximation.” Brocot tells us that this section was “‘spe- 
cially intended for English clock-makers.” Compare RMT 87. 

R. C. A. 


CORRIGENDA ET ADDENDUM 


P. 6, 54, for circular and cosine, read circular sine and cosine. 

P. 19, 1. 15 from bottom, for MARTHA JAEGER, read MARTHA CLARKE. P. 33, |. 11 and p. 100, 
for THEODERICUs, read THEODoricus. P. 40, 1. 3, for 639000, read 63900. P. 42, 1. 6, for 
1848, read 1858. P. 44, footnote 7, 1. 1, for g, , “, read g, ‘, “. P. 45, 1. 26, for 1,0000, 
read 1.0000; for 1.000, read 1.0000; for 25000, read 2.5000. 

P. 46, 1. 10, for 0:00100, read 0.00100. P. 50, for e-#, read e~#*; for dx, read dx"; for h(x), 
read H(x); for and H(0.3), read and H(0.2). P. 51, 1. 7, for Levinson, read LEVENSON, 
P. 55, 1. 7, for 6.974, read 6.474. P. 61 MAC, 1. 11, for specializing, read specialize. P. 
70, 1909, 1. 2, read 100(100)1000; 1912 for specimal, read specimen. P, 71, 1914, 1. 3 read 
C,’ = dC,(x)/dx. P. 72, 34, for A. E. KENNELLY, read P. L. ALGER. 

P. 74, 54, 1. 3, read J3(x)/Is(x). P. 75, 1. 5 from bottom, for .6(.2)1, read .6(.02)1. P. 76, 
1. 11, for (.005)1(.001)-9, read (.005).1(.001).9. P. 84, (7.c.s.), read . . . huit Décimales 
des Nombres entiers de 1 2 120000 et des Sinus et Tangentes de dix Secondes en dix 
Secondes d’Arc . . .; (11 t.s.), 1. 3, for 1%, read 10°. P. 75, for KENNELLY, read ALGER. 
P. 91, Table III, 1. 5, for —7.0*, read —70*. 

P. 94, 1. 5, from bottom for 0(.001).99, read 0(.0001).99. P. 95, 1. 5, for 5D, read to degrees and 
minutes; |. 13, for 4S, read 4D; 1. 20, read Umkehrfunktionen. P. 96, for 5(X), read 5(Z). 
P. 99, 1. 13, for 0(0.01)3, read 0(0.1)3. 

P. 69, 1879, add to the described contents of Glaisher’s article as follows: On p. 53 are given 

the zeros of P,(x), for m = [2(1)7; 15D]. 
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